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ABSTRACT 


The  object  of  this  thesis  is  to  develop  techniques  for 
the  iterative  solution  of  problems  on  the  analog  computer. 

An  essential  requirement  for  iterative  computation  is  storage. 
For  this  purpose  an  iterative  control  unit  was  built  contain¬ 
ing  the  necessary  logic  and  solid-state  switching.  An  analysis 
is  made  of  the  various  types  of  memories  that  may  be  imple¬ 
mented  using  solid-state  switching  on  the  analog  computer. 

The  application  of  iterative  techniques  to  the  solution  of  a 
two-point  boundary  value  problem  is  given  for  a  linear  second 
order  differential  equation  for  illustrative  purposes.  The 
generalization  to  more  complicated,  nonlinear  problems  is 
immediate  since  the  iteration  scheme  remains  essentially 
unchanged.  In  the  serial  solution  of  partial  differential 
equations,  the  value  of  the  function  from  the  previous  inter¬ 
val  of  time  is  required.  A  technique  is  developed  where  this 
storage  may  be  accomplished  with  only  moderate  equipment 
requirements.  This  makes  the  serial  solution  of  partial 
differential  equations  using  an  analog  computer  more  realistic 
since  function  storage  normally  imposes  excessive  equipment 
requirements.  Using  this  function  storage,  the  serial  solution 
to  a  parabolic  partial  differential  equation  is  obtained. 
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INTRODUCTION 


1  .0 

Iterative  computation  represents  one  of  the  most 
effective  implementations  of  an  analog  computer.  In  normal 
analog  computing,  many  trial  solutions  are  made  before  the 
final  set  of  parameters  are  chosen.  This  trial  technique 
may  be  accelerated  if  the  solutions  to  the  programmed 
equations  are  obtained  in  repetitive  operation.  Since  the 
solutions  are  displayed  visually,  the  correct  selection 
of  variables  is  more  direct.  The  repetitive  feature  how¬ 
ever  itself  only  provides  a  means  for  the  visual  display  of 
results.  In  iterative  computation,  the  solutions  are  also 
obtained  in  repetitive  operation.  However  the  results  obtained 
in  each  operate  cycle  are  used  to  change  some  condition  such 
as  an  initial  condition  or  a  parameter  so  as  to  make  the 
next  solution  closer  to  that  desired.  That  is,  the  procedure 
of  adjusting  various  parameters  for  the  solution  to  meet 
specified  conditions  has  been  automated. 

The  reason  iterative  computation  may  be  referred  to  as 
the  most  effective  implementation  of  an  analog  computer  may 
be  exemplified  by  the  consideration  of  a  two-point  boundary 
value  problem.  The  normal  procedure  would  be  to  perform 
many  runs,  changing  the  required  initial  condition  until 
such  time  that  the  boundary  conditions  were  satisfied.  This 
procedure  is  both  lengthy  and  boring.  This  same  problem 
using  iterative  techniques  can  be  solved  automatically  in 


about  60  milliseconds. 
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Iterative  techniques  may  be  applied  effectively  to 
many  problems  on  the  analog  computer.  Some  of  the 
applications  are  listed  below. 

(a)  Boundary  Value  Problems.  In  this  application  an 
iteration  is  performed  to  find  a  remaining  initial 
condition.  Convergence  to  the  required  result  is 
both  fast  and  reliable. 

(b)  Optimization  Problems.  In  this  case,  some  performance 
criterion  such  as  the  integral  of  the  error  squared 

of  a  system  may  be  minimized  by  varying  one  or  more 
parameters.  The  iterative  technique  also  faciliates 
the  application  of  Pontryagin’s  Maximum  Principle 
to  optimization  problems  on  the  analog  computer. 

(c)  Solution  of  Partial  Differential  Equations.  The 
iterative  technique  may  be  applied  to  obtain  the  serial 
solution  of  a  partial  differential  equation.  Equipment 
requirements  become  more  severe  in  this  application 
because  of  the  necessity  for  function  storage. 

(dj  Parameter  Sweep.  In  this  application  iterative 

techniques  may  be  aoplied  to  obtain  for  instance,  a 
plot  of  the  peak  overshoot  of  a  system  as  a  function 
of  f .  The  proper  selection  of  f  is  then  immediate. 

This  same  technique  may  also  be  used  to  solve 
optimization  problems  where  a  performance  functional 
is  to  be  minimized  with  respect  to  one  parameter,  say 
<x.  This  scheme  can  then  be  used  to  obtain  a  plot  of 
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this  functional  as  a  function  of  ex-  on  the  x-y  plotter.  The 
optimum  choice  of  is  then  immediate. 

ie)  Curve  Fitting.  Iterative  techniques  may  be  used  to 

extrapolate  a  curve  through  a  given  set  of  points  so 
so  as  to  satisfy  some  error  criterion.  This  is  similar 
to  the  case  where  a  nerformance  functional  is  being 
minimized  in  a  optimization  problem. 

(f)  Solution  of  Eigenvalues.  Iterative  techniques  may 
be  used  on  the  analog  computer  to  find  the  first  set  of 
eigenvalues  in  boundary  value  problems  in  mathematical 
physics.  The  approach  is  similar  to  that  for  the  two-point 
boundary  value  problem,  only  in  this  case  the  eigenvalue 
would  be  the  parameter  rather  than  the  initial  condition. 

1.1  CONTRIBUTIONS  OF  THIS  THESIS 

In  this  thesis  the  application  of  iterative  techniques 
to  an  analog  computer  are  investigated  using  the  Face  231 -R 
general  purpose  computer.  An  iterative  control  unit  was  built 
containing  the  necessary  logic  and  solid-state  switching 
required  for  iterative  computation.  The  implementation  of 
solid-state  switching  enabled  iterative  computations  to  be 
performed  at  higher  speeds  with  improved  accuracy  as  compared 
to  normal  relay  switching. 

The  first  section  deals  with  the  development  of  solid- 
state  switching  techniques  to  perform  storage  using  the 


analog  computer.  An  analysis  of  the  various  types  of 
memories  and  their  respective  performance  is  given. 

The  iterative  solution  of  a  two-point  boundary-value 
problem  was  applied  to  a  second  order  ordinary  differential 
equation  for  illustrative  purposes.  The  generalization  to 
more  complicated  problems,  for  instance,  nonlinear  problems 
is  immediate  since  the  iteration  scheme  remains  essentially 
unchanged.  A  very  significant  application  of  the  iterative 
solution  of  a  two-point  boundary  value  problem  is  in  the 
application  of  Pontryagin's  Maximum  Principle  to  optimization 
problems  on  the  analog  computer. 

The  major  objective  of  this  thesis  was  to  perform  the 
serial  solution  of  a  partial  differential  equation,  using 
solid-state  switching.  The  serial  solution  of  a  partial 
differential  equation  however  imposes  severe  equipment 
requirements  because  of  the  need  for  function  storage.  A 
method  was  developed  where  the  equipment  requirement  was 
reduced  from  seventy  amplifiers  to  just  twenty-five  amplifiers 
for  an  automatic  function  generator  using  data  stored  for 
nine  equally  spaced  intervals  of  the  independent  variable. 

With  the  development  of  mode  control  using  field  effect 
transistors,  this  requirement  may  be  reduced  to  just  fifteen 
amplifiers  for  parabolic  curve  fitting  using  any  number  of 
stored  points.  The  program  of  extrapolating  the  function 
through  the  stored  points  is  referred  to  as  Newton's 
Iterative  Function  Generator.  Since  this  function  generator 


requires  but  a  moderate  number  of  amplifiers  for  function 
storage,  tnis  scheme  can  be  applied  to  smaller  computer 
installations,  for  instance  the  race  TR-lj.8  computer.  An 
illustration  of  the  apnlication  of  the  serial  method  to  the 
solution  of  a  parabolic  partial  differential  equation  is 
given  complete  with  an  analysis  of  the  errors  involved. 

1.2  A  SURVEY  OP  THE  CHAPTERS 

Chapter  2  Tne  technique  for  implementing  storage  on 

the  analog  computer  using  solid-state  switching  is 
developed.  A  description  of  the  various  types  of 
memories  that  are  used  in  iterative  analog  computation 
is  given. 

Chapter  3  The  iterative  solution  of  a  two-point  boundary 

value  problem  is  given.  A  second  order  ordinary 
differential  equation  is  used  for  illustrative  purposes. 
Chapter  A  technique  is  developed  for  the  storage  of 

a  function  with  only  moderate  equipment  requirements. 

A  program  is  given  for  an  automatic  i unction  generator 
using  nine  stored  points. 

Chapter  5  Using  the  function  storage  developed  in 

Chapter  I4.,  the  serial  solution  to  a  parabolic  partial 
differential  equation  is  given.  An  analysis  of  the 
incountered  in  this  method  is  also  given. 
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2.U  STORAGE  ON  THE  ANALOG  COMPUTER 


The  basic  requirement  in  the  application  of  iterative 
techniques  to  problems  on  the  analog  computer  is  that  of 
storage.  There  are  a  variety  of  ways  in  which  this  may  be 
accomplished.  One  is  the  use  of  a  Memory  Wheel  as  described 
by  Kozak  (reference  lOj.  However  this  technique  is 
mechanically  cumbersome  and  syncroniza tion  with  a  computer 
in  repetitive  operation  becomes  a  serious  problem. 

A  second  technique  that  can  be  employed  is  that  of 
Hybrid  Computation.  A  small  digital  computer  could  perform 
necessary  storage  as  well  as  some  sub-routines.  Such  a 
facility  however  was  not  available. 

The  third  technique  and  the  one  used  in  this 
investigation  is  storing  a  charge  on  a  capacitor  across  a 
switched  amplifier.  This  technique  lends  it  self  very  well 
to  cases  where  the  storage  requirements  are  moderate. 

Another  favorable  aspect  of  this  scheme  is  that  only  a 
modest  amount  of  extra  equipment  is  required. 

2.1  STORAGE  ON  AN  OPERATIONAL  AMPLIFIER 

The  basic  mechanism  employed  in  storage  using  an 
operational  amplifier  is  given  in  Figure  2-1.  When  switch 

A  is  closed, 

=  J . . (2-1  ) 

ein  1  +rD 


where  'f—  RC 
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ein 


eout 


FIGURE  2-1.  SCHEMATIC  FOR  STORAGE  AMPLIFIER 


If  r  is  made  small  with  respect  to  the  time  constant  in 
the  system  being  solved,  we  may  approximate  the  expression 
(2-1 )  by, 


When  switch  A  is  open,  the  output  voltage  will  remain 
constant.  That  is  the  amplifier  will  now  hold  or  store 
the  value  at  the  switching  instant. 

Switch  A  is  a  diode  bridge  driven  by  a  set  of 
complementary  flip-flops,  which  are  triggered  by  a  Schmitt 
Trigger  comparator.  A  schematic  for  the  switching 
arrangement  is  given  in  Figure  2-2. 

In  the  amplifiers  used  for  storage, 

R  =  1 00K  ohms 
C  =  .001  ufd. 

Therefore  fT/  —  100  x  10^  xl  xlo  '  =  1  xlo  seconds. 

The  time  required  for  an  amplifier  to  track  a  voltage  to 


within  1  percent  error  is; 
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t  -  U.  61  -  0 . 1.j.6  milliseconds 
1 

T^is  time  constant  can  be  made  smaller  by  using  a  .0U01  ufd. 
capacitor.  However  a  comprimise  must  be  made  between  the 
ability  for  a  memory  to  track  a  voltage  accurately  and  the 
ability  for  this  memory  to  hold  or  store  a  value. 


FIGURE  2-2.  SWITCHING  FOR  STORAGE  AMPLIFIER 


2.2  FORWARD  AND  REVERIE  MEMORIES 

A  forward  memory  (FMJ  is  one  that  tracks  when  the 
computer  is  in  the  operate  mode  and  stores  or  holds  its 
value  when  the  computer  is  in  the  hold  or  reset  mode. 

A  reverse  memory  (RMj  is  the  complement  of  a  forward 
memory.  It  will  track  when  the  computer  is  in  the  reset 
mode  and  will  store  when  the  computer  is  in  the  operate 
mode.  The  symbols  for  each  memory  are  given  in  Figure  2-3. 
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FIGURE]  2-3.  FORWARD  AND  REVERSE  MEMORY  SYMBOLS 


A  graphical  explanation  of  the  performance  of  each 
memory  is  given  in  Figures  2-li  and  2-5> . 


FM  input 
FM  output 


FIGURE  2-k-  OPERATION  OF  A  FORWARD  MEMORY 
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RM  input  - 

RM  output  - 

FIGURE  2-5.  OPERATION  OF  A  REVERSE  MEMORY 

2.3  ELECTRONIC  COMPARATOR 

The  purpose  of  the  electronic  comparator  in  storage 
techniaues  on  the  analog  computer  is  to  give  flexibility 
to  the  program  by  making  it  possible  to  switch  the  mode 

of  a  memory  at  a  particular  instant  during  the  operate  cycle. 
A  schematic  of  the  electronic  comparator  is  given  in 

Figure  2-6.  The  symbols  N  and  S  refer  to  the  outputs  of 
the  comparator  in  the  normal  and  switched  state  respectively. 


input- 


Schmitt 

n-p-n 

Trigger 

nip- 
flop  | 

+12N 

- .  63 


i 

i 


A 


FIGURE  2-6.  ELECTRONIC  COMPARATOR 
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The  symbols  that  will  be  used  for  designating  electronic 
comparators  are  given  in  Figures  2-7  and  2-8. 


FIGURE  2-7.  ELECTRONIC  COMPARATOR  SYMBOL 


Figure  2-7  is  the  symbol  for  the  electronic  comparator 
described  above.  At  the  instant  the  input  voltage  exceeds 
the  negative  reference,  the  state  of  the  comparator  will 
change  and  the  new  output  voltages  will  be  as  indicated 
for  the  switched  state  in  Figure  2-6. 

When  the  number  of  input  voltages  to  be  compared  is 
greater  than  two,  a  unity  gain  amplifier  is  used  as  a  summer 
for  the  comparator.  The  symbol  that  will  be  used  in  this 
case  is  shown  in  Figure  2-8. 


FIGURE  2-8.  ELECTRONIC  COMPARATOR  SYMBOL 


In  Figure  2-8,  B  designates  the  summing  amplifier 
while  M-j  refers  to  the  number  of  the  electronic  comparator. 
A  latch  voltage  (El)  is  sometimes  used  which  is  larger 
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than  the  signal  and  reference  voltages.  The  state  of  the 
comparator  will  then  be  determined  only  by  the  sign  of  the 
latch  voltage.  When  the  signal  voltage  is  less  than  the 
reference  voltage  (and  SL  =  0;,  the  output  voltage  from  the 
comparator  will  be  identical  to  that  shown  in  Figure  2-6 
for  the  normal  state.  At  the  instant  the  signal  voltage 
exceeds  tne  positive  reference  voltage,  the  output  voltages 
A  and  A  will  assume  the  values  given  in  Figure  2-6  for  the 
switched  state. 

2.k  X-MEMORIES  AND  0-MEMORIES 

X-Memories  and  0-Memories  perform  essentially  the  same 
function  as  the  forward  and  reverse  memories,  however  they 
are  comparator  controlled.  An  extra  degree  of  flexibility 
is  added  when  using  this  type  of  memory  for  storage  since 
at  any  particular  instant  in  the  operate  cycle  a  value  may 
be  stored. 

The  function  of  an  x-memory  is  related  to  a  forward 
memory  in  that,  for  the  first  portion  of  the  operate  cycle 
this  memory  will  track  its  input  voltage.  At  some  particular 
instant  during  the  onerate  cycle  the  mode  of  this  memory 
will  be  changed  by  an  electronic  comparator,  and  this  memory 
will  now  hold  its  last  value.  The  symbol  for  an  x-memory 
in  given  in  Figure  2-9. 
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The  x-memory  in  Figure  2-9  will  track  the  input  voltage 
until  x  =  L,  at  which  time  its  mode  will  change.  The  memory 
will  now  store  or  hold  the  value  f(L)  for  the  remainder 
of  the  operate  cycle.  Mode  control  is  furnished  by  comparator 
in  the  following  fashion,,  When  the  input  signal  is  less 
than  the  negative  reference,  the  comparator  is  in  its  normal 
state.  In  this  normal  state,  output  A  will  cause  the  diode 
bridge  used  in  the  x-memory  to  conduct.  This  will  cause 
the  memory  to  track  its  input  voltage.  At  the  instant 
x  =  L,  the  state  of  the  comparator  will  change  and  output 
A  will  produce  voltages  so  as  to  make  the  diode  bridge  in 
the  mempry  non-conducting.  This  in  turn  will  cause  the 
the  x-memory  to  store  the  value  of  the  input  function  at 
x  =  L.  Figure  2-10  displays  the  operation  of  an  x-memory. 

An  0 -Memory  is  the  complement  of  an  x-memory  in  the 
same  sense  the  reverse  memory  was  the  complement  of  a  forward 
memory.  The  symbol  for  an  o-memory  is  given  in  Figure  2-11. 


-  1*1.  - 


X-M  input  - 

X-M  output  - 

FIGURE  2-10.  OPERATION  OF  X-MEMORY 


+  input  signal 
-100 


eout 


FIGURE  2-11.  O-MEMORY 


When  comparator  M^  is  in  the  normal  state,  voltage 
A  will  cause  the  diode  bridge  in  the  o-memory  to  non-conduct. 
In  this  mode  the  o-memory  will  hold  or  store c  At  the  instant 
the  input  signal  is  equal  to  L,  the  comparator  will  switch 
its  state  and  voltage  A  will  in  turn  cause  the  bridge 
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to  conduct.  The  o-memory  will  now  track  its  input  voltage. 
Figure  2-12  displays  the  operation  of  an  o-memory. 


Computer  Mode 

0-M  input  - 

0-M  output  - 

FIGURE  2-12.  OPERATION  OF  AN  O-MEMORY 

2.5  RATCHET  MEMORY 

By  cascading  an  x-memory  and  an  o-memory  a  value  can 
be  remembered  from  a  previous  cycle.  This  technique  is 
the  key  to  iterative  computation  on  the  analog  computer. 

An  example  where  a  ratchet  memory  can  be  used  is  in  the 
solution  of  an  optimization  problem  where  some  parameter  or 
parameters  are  being  varied  in  order  to  minimize  some 
performance  functional.  The  value  of  the  performance 
functional  from  the  previous  iteration  is  required  to  make 
corrections  in  these  parameters  in  order  to  eventually  arrive 
at  the  optimum  choice.  A  schematic  of  a  ratchet  memory  is 
given  in  Figure  2-13. 
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FIGURE  2-13.  RATCHET  MEMORY 

During  the  first  portion  of  the  operate  cycle,  the 
x-memory  will  track  the  input  voltage  f(t)  while  the  o-memory 
remains  in  the  hold  mode.  At  the  instant  x  =  L,  the  x-memory 
will  store  its  last  value,  that  is  -f(L),  while  the  o-memory 
will  now  read  in,  or  track  the  voltage  held  by  the  x-mernory. 

At  the  completion  of  the  operate  cycle,  the  x-memory  will 
again  track  its  input  voltage  while  the  o-memory  will  hold 
or  store  the  value  f(Lj.  At  the  beginning  of  the  next  operate 
cycle  the  o-memory  will  still  be  storing  the  value  f(L).  In 
this  fashion  a  value  has  been  remembered  from  a  previous 
computation.  A  number  of  ratchet  memories  may  be  cascaded 
so  that  in  the  n'th  computation,  values  of  some  functional 
or  parameter  from  the  n-1 ,  n-2,  . . .  computation  may  be 
available.  The  schematic  for  a  cascaded  ratchet  memory  is 
given  in  Figure  2-1 1|.  The  figure  is  drawn  for  the  n'th 
computation  with  values  of  f(L)  for  the  last  two  computations 
being  made  available.  The  generalization  to  the  storage  of 
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this  value  from  more  computations,  for  instance  n- 3 ,  n-lp,  .  .  , 
is  immediate. 


FIGURE  2-1 1^.  CASCADED  RATCHET  MEMORY 

2.6  ACCURACY  OF  MEMORIES 

A  memory  using  a  .001  ufd.  capacitor  for  storage  will 
track  a  fixed  voltage  to  within  *0^  volts.  A  further 
improvement  could  be  made  by  balancing  the  diode  bridges 
to  closer  tolerances.  In  the  hold  state,  drift  was  normally 
less  than  0.1  volts  in  ten  seconds.  This  could  be  further 
improved  by  placing  an  external  bias  to  compensate  for  any 
bridge  unbalance.  In  this  manner  drift  in  the  hold  state 
is  less  than  0.01  volts  in  ten  seconds.  When  using  a  larger 
capacitor  for  storage,  error  in  the  track  and  hold  modes 
was  imperceptible. 


3.0  SOLUTION  OP  BOUNDARY  VALUE  PROBLEMS  IN  ORDINARY 
DIFFERENTIAL  EQUATIONS 

The  technique  of  storage  on  the  analog  computer  facili¬ 
tates  the  automatic  solution  of  a  boundary  value  problem  in 
ordinary  differential  equations.  Only  one  ratchet  memory 
is  required  when  one  boundary  value  is  being  sought,  A 
typical  problem  is  a  trajectory  problem  (reference  6j,  where 
the  firing  angle  is  to  be  determined  for  a  missile  which  is 
to  strike  a  target  a  distance  A  away.  One  method  for  solving 
this  problem  is  to  program  the  applicable  equations  on  the 
analog  computer,  then  run  the  solution  in  repetitive  operation. 
By  manual  adjustment,  the  required  firing  angle  may  be  found. 
However,  by  the  addition  of  a  simple  ratchet  memory  this 
process  may  be  automated. 

In  the  solution  for  an  unknown  boundary  value  in  linear 
ordinary  differential  equations,  the  advantage  of  iterative 
computation  is  not  evident  with  the  exception  of  illustrative 
purposes.  If  the  problem  Is  made  nonlinear,  for  instance,  if 
air  friction  is  taken  into  consideration  in  the  missile 
problem,  then  the  advantage  of  iterative  computation  becomes 
more  evident. 

In  the  serial  solution  of  a  partial  differential  equation 
on  an  analog  computer,  for  each  instant  of  time  a  two-point 
boundary  value  problem  in  ordinary  differential  equations 
must  be  solved  (referred  to  as  a  nested  calculation).  The 
iteration  scheme  required  in  the  solution  of  a  partial 
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differential  equations  is  practically  identical  to  that  in 
in  the  solution  of  a  boundary  value  problem  in  ordinary 
differential  equations.  An  illustration  of  the  application 
of  iterative  techniques  to  a  boundary  value  problem  in  a 
second  order  differential  equation  will  be  given  in  Section 
3-1  . 


3.1  BOUNDARY  VALUE  PROBLEM  IN  ORDINARY  DIFFERENTIAL 
EQUATIONS 


The  technique  of  applying  storage  techniques  to  obtain 
solutions  to  boundary  value  problems  will  be  illustrated 
by  considering  a  simple  second  order  differential  equation 
(reference  11).  The  problem  that  will  be  considered  is 
given  below. 

Given  the  O.D.E. 

d9y  +  y  =  0  . (3-1  ) 

dx^ 

Find  d£  such  that  the  following  boundary  conditions 

dx  x=0 

are  satisfied. 


y (0 )  =  0 

y  ( L )  =  A  0<L<rT  A  <  1 

The  theoretical  solution  may  readily  be  obtained.  This  is 
equal  to, 


y(x) 
y  (L) 


A  sin  x 


therefore 


A 


A  sin  L 

y (L)  . . 
sin  L 


(3-2) 
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3.1.1  COMPUTER  PROGRAM 


The  necessary  scaling  will  first  be  performed.  The 
normalized  variables  that  will  be  used  as  as  follows: 


iz  =  1 _ _  rd£l 

dx  100  dx 


y 


x 


1  rx  ~| 
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The  computer  program  is  given  in  Figure  3-1 .  Scaling 
in  the  x  variable  is  for  10  millisecond  repetitive  operation. 
Integrators  1 u  and  11,  and  amplifier  12  comprise  the  program 
for  solving  equation  (3-1).  This  equation  is  solved  in 
repetitive  operation  on  the  computer.  For  the  first  portion 
of  the  operate  cycle,  that  is  when  x<  L,  amplifier  17  will 
track  the  voltage  A  -  y(x).  At  the  instant  x  =  L,  comparator 
1  will  switch  state,  changing  the  mode  of  x-memory  17  and 
o-memory  18  so  the  x-memory  17  will  hold  the  value  A  -  y(L) . 
0-Memory  18  will  now  track  this  voltage.  At  the  end  of  tne 
operate  cycle  the  mode  of  the  memories  will  again  change 
and  o-memory  18  will  now  store  tne  value  A  -  y(L).  But 
A  -  y(L)  serves  as  a  practical  value  for  the  error  in  a  given 
iteration.  Integrator  13  will  now  integrate  at  a  rate 
proportional  to  this  error  during  tne  reset  cycle.  Amplifier 
14  will  now  provide  a  better  approximation  for  the  required 
initial  condition.  This  iteration  process  will  continue  until 
tne  required  initial  condition  has  been  obtained. 


- 
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FIGURE  3-1 .  COMPUTER  PROGRAM 
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During  the  operate  cycle,  integrator  13  will  remain  in 
the  hold  mode.  This  mode  control  is  performed  by  amplifiers 
07,  08,  and  diode  switch  B.  A  detailed  description  for 
external  mode  control  for  integrators  is  given  in  Appendix  P. 
Potentiometer  1)4.  provides  an  initial  guess  for  the  required 
initial  condition. 

A  photograph  showing  the  successive  iterations  to  the 
solution  is  given  in  Figure  5-12 a.  In  this  particular  case 
L  =  ^ / 2  and  y(L)  was  equal  to  0.5.  The  results  obtained  gave 
A  =  .502  .  This  is  an  error  of  less  than  one  half  of  one 
percent,  making  the  iterative  technique  of  computation 
valuable  in  the  solution  of  two-point  boundary  value  problems. 

3.2  CONCLUSIONS 

The  application  of  iterative  techniques  on  the  analog 
computer  proves  to  be  a  very  effective  method  of  solving 
two-point  boundary  value  problems  in  ordinary  differential 
equations.  If  the  equations  are  further  complicated  by 
nonll neari ties ,  the  only  required  addition  may  be  a 
multiplier  and  a  diode  function  generator.  The  iteration 
scheme  however  remains  essentially  identical  to  that  for  the 
linear  case. 

This  technique  of  solving  two-point  boundary  value 
problems  may  be  applied  to  optimization  problems  employing 
Pontryagin's  Maximum  Principle.  The  possibility  of  using 
field  effect  transistors  with  Nexus  Operation  amplifiers 
for  repetitve  computations  at  a  frequency  of  30  kilocycles 
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per  second  has  been  investigated  in  the  Electrical  Engineering 
Department  of  the  University  of  Alberta  with  very  promising 
results  (reference  -Ilf).  This  would  imply  that  iterative 
computations  at  a  frequency  of  at  least  10  kilocycles  could 
easily  be  performed.  This  would  present  an  opportunity  for 
the  effective  application  of  Pontryagin* s  Maximum  Principle 
to  "in-line"  optimization  problems  on  the  analog  computer. 


U*0  FUNCTION  STORAGE  ON  THE  ANALOG  COMPUTER 


Function  storage  is  required  when  the  serial  solution 
of  a  partial  differential  eauation  is  obtained  using  an 
analog  computer.  In  this  particular  case,  the  value  of  the 
solution  from  the  previous  interval  of  time  is  required 
for  obtaining  the  solution  for  the  following  interval  of 
time.  The  need  for  curve  storage  occurs  in  other  situations 
such  as  in  optimization  problems  and  in  automatic  curve 
fitting  on  the  analog  computer. 

Curve  storage  on  the  analog  computer  imposes  stringent 
equipment  requirements.  In  the  case  of  the  serial  solution 
of  a  parabolic  partial  differential  equation  approximately 
eighty  percent  of  the  amplifiers  available  in  the  computer 
are  used  for  implementing  curve  storage  and  read  out. 

A  block  diagram  for  the  general  technique  for  function 
storage  on  the  analog  computer  is  given  in  Figure  I4.-I  . 


FIGURE  1±-1  .  BLOCK  DIAGRAM  FOR  FUNCTION  STORAGE. 


The  sampler  consists  of  a  series  of  comparator  controlled 
x-memor ies  „  One  of  these  x-memories  is  given  in  Figure  I+-2. 
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FIGURE  \\-2. .  X-MEMORY  IN  SAMPLER 


X-Memory  M  will  track  the  input  voltage  f(t)  until  such 
time  that  t  =  tm,  when  the  comparator  will  switch  the  mode 
of  M  and  the  x-memory  will  now  hold  or  store  the  value  of 
f(tm).  In  a  similar  fashion  other  x-memories  will  store 
the  values  f(t0),  f  ( t  >|  ) ,  f(t2),  ,*...f(tn).  In  all  further 
analysis  the  function  will  be  considered  to  be  sampled  at 

t 

equal  intervals  of  time. 

The  automatic  function  generator  will  consist  of  a 
program  on  the  analog  computer  such  that  using  the  stored 
values,  f(t0),  f  ( t  -j  )  ,  .....f(tn),  the  continuous  function 
f(t)  can  be  regenerated  when  the  program  is  driven  by  a 
voltage  proportional  to  time.  The  automatic  function 
generator  will  employ  the  Newton  Forward  Difference 
Interpolation  formula  for  regeneration  of  the  required 
function.  The  Newton  Forward  Difference  Formula  is; 

f(t)  =  a  +  bt  +  ct(t-l)  +  dt(t-1)(t-2)  +  ... 

2!  3* 

Coefficients  a,  b,  c,  and  d  are  obtained  in  Table  I4.-I 


( reference  1  )  . 
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0 

1 

2 

3 


f  (t) 

Ko 

Kl 

k2 

Ko 


ZK 


K1  -  K0 
K2  -  K>, 

-  k2 


Z 


3 


K2  -  2K>|  +  Kq 

Ko  -  3K2  +  3K>] 
K3  -  2K2  +  3L,  -  Ko' 


TABLE  4-1.  COEFFICIENTS  FOR  NEWTON’S 

FORWARD  DIFFERENCE  FORMULA. 


From  Table  4“ 1  ,  it  can  be  deduced  that, 

a  =  K0 

b  =  K,  -  K0 

c  =  K2  -  2K/|  +  K0 

d  =  K3  -  3K2  +  3K>,  -  K0 

This  procedure  for  finding  the  coefficients  can  easily  be 
generalized  to  the  case  where  more  points  are  being  used 
(See  Appendix  G-). 


4.1  PROGRAMMING  FOR  AN  AUTOMATIC  FUNCTION  GENERATOR 

There  are  three  schemes  which  can  be  used  in  programming 
an  automatic  function  generator,  A  brief  discription  of  each 
will  be  given. 


4.1.1  METHOD  1.  (DIRECT  SIMULATION) 

The  scheme  shown  in  Figure  4-3  is  that  for  three  stored 
points  (reference  $).  The  generalization  to  more  points 


27  - 


is  immediate. 


.01 


Using  the  above  scheme  for  nine  points  would  require 
sixty  amplifiers  and  fifteen  multipliers.  As  is  evident, 
this  scheme  could  be  used  successfully  for  up  to  about  four 
points.  However  if  better  accuracy  is  required  and  hence 
more  points  must  be  used  in  the  curve  extrapolation,  the 
excessive  requirement  for  multipliers  makes  this  scheme 
impractical . 

J+.1.2  METHOD  2.  ( INDIRECT  SIMULATION ) 

A  block  diagram  of  the  scheme  employed  is  given  in 
Figure  ip— ip  (reference  7). 
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FIGURE  4-4.  AUTOMATIC  FUNCTION  GENERATOR  (INDIRECT  METHOD). 

Using  the  table  drawn  up  in  Appendix  G,  an  algebraic 
calculation  can  be  performed  on  the  computer  to  find  the 
coefficients  a,  b,  c,  .  .  .  .n.  For  example,  b  =  K-j  -  K0  . 

If  the  function  and  its  derivatives  are  considered  at  zero 
we  have; 

f(t)  =  a  +  bt  +  jrt(t-l)  +  ... 

2! 

f(0)  =  a 

f'(t)  =  b  +  ct 

f'(0)  =  b 

f  \  t )  =  c 

f(0)  =  c 

After  all  the  algebraic  calculations  for  the  values  a,  b,  c 
are  completed,  a  second  algebraic  calculation  can  be  performed 
to  give  the  value  of  the  function  as  well  as  all  of  its 
derivatives  at  t  =  0.  Examining  the  equation  more  carefully 
shows  that  the  highest  non-zero  derivative  of  the  polynomial 
gives  the  derivative  of  the  function  throughout  the  entire 
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period  of  interest.  In  the  above  case,  f(t)  =  c.  Since 

a 

l(t)  and  all  initial  values  of  the  derivatives  of  all  lower 
orders  are  known,  successive  integration  may  be  applied  to 
obtain  the  desired  function. 

li.1.3  METHOD  3.  (NEWTON'S  ITERATIVE  FUNCTION  GENERATOR) 

The  technique  employed  here  is  to  use  two  sets  of  Newton's 
Function  Generators,  each  programmed  for  n  points  using  the 
technique  described  in  Method  2.  Each  is  switched  in  or  out 
of  operation  by  a  comparator,  hence  giving  a  total  possible 
number  of  2n  stored  points.  A  block  diagram  is  given  in 
Figure  4"  5 . 

This  technique  can  easily  be  generalized  to  provide 
a  means  for  extrapolating  stored  data  from  a  digital  memory. 

The  stored  function  could  then  be  read  out  in  the  form  of  a 
polynomial  rather  than  in  box-car  form.  This  is  depicted 
graphically  in  Figure  4“5a  . 


FIGURE  4-5.  BLOCK  DIAGRAM  FOR  NEWTON'S  ITERATIVE 


FUNCTION  GENERATOR. 
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Solved  by  A  Solved  by  B  Solved  by  A 
(B  in  I.C.)  (A  in  I.C.)  (B  in  I.C.) 

FIGURE  4-. 3a.  NEWTON* S  ITERATIVE  FUNCTION  GENERATOR 
FOR  DIGITAL  TO  ANALOG  CONVERTER 

A  brief  description  of  the  operation  of  Newton’s 
Iterative  Function  Generator  as  a  digital  to  analog  converter 
is  given  with  reference  to  Figure  1+-3 .  When  program  A  is  in 
the  operate  mode,  B  is  in  the  initial  condition  mode  with  its 
output  switch  openc  At  the  end  of  operate  cycle  of  program 
A,  the  output  switch  of  A  will  open,  that  of  B  will  close  and 
3  will  now  be  in  the  operate  mode.  This  sequence  is  then 
repeated  to  extrapolate  the  entire  curve. 

NEWTON’S  PROGRAM  FOR  NINE  POINTS 

To  program  an  automatic  function  generator  for  nine 
points  requires  twenty-five  amplifiers  using  Newton’s  Iterative 
Method  as  compared  to  sixty  amplifiers  and  fifteen  multipliers 
for  Method  1,  and  seventy  amplifiers  for  Method  2.  The 
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advantage  of  employing  this  technique  is  evident,, 

A  memory  programmed  for  9  points  will  be  used  in  the 
serial  solution  of  a  parabolic  partial  differential  equation. 
Two  automatic  function  generators,  each  of  five  points  are 
programmed  and  are  switched  in  the  final  integrator,  that  is, 
the  integrator  generating  -f(tj.  This  is  done  because  an 
attempt  to  perform  switching  in  the  summing  junction  of  an 
amplifier  tends  to  produce  a  ringing  effect.  A  block 
diagram  for  an  automatic  function  generator  using  9  stored 
points  is  given  in  Figure  J4.-6. 


FIGURE  4-6.  FUNCTION  GENERATOR  BLOCK  DIAGRAM  (9  POINTS). 

The  method  of  operation  of  the  program  in  Figure  4“ 6  is 
similar  to  that  described  for  the  Newton’s  Iterative  Function 
Generator.  However  the  initial  conditions  are  read  in 
simultaneously  to  both  programs  when  the  computer  is  in  the 
reset  mode.  Program  A  then  generates  the  function  for  the 
first  four  segments  during  which  time  program  B  is  forced  to 
remain  in  the  initial  condition  mode  (with  switch  B  open). 
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At  x  =  L/2,  switch  A  is  opened  and  switch  B  is  closed. 

Program  B  now  generates  the  remaining  half  of  the  function. 

No  extra  initial  condition  is  required  for  B  since  the  desired 
value  or  initial  condition  for  program  B  is  already  existing 
on  integrator  1  when  A  completes  its  operate  cycle.  Program 
A  is  forced  to  remain  in  the  hold  mode.  Otherwise  program  A 
would  continue  the  generation  of  some  extension  of  the  first 
half  of  the  desired  curve  and  eventually  saturate  some  ampli¬ 
fiers.  The  technique  employed  for  having  B  remain  in  the 
I.C.  mode  for  the  first  half  of  the  operate  cycle  and  to 
have  A  remain  in  the  hold  mode  for  the  remainder  of  the 
operate  cycle  is  described  using  Figure  it-7. 


FIGURE  li-7.  MODE  CONTROL  FOR  FUNCTION  GENERATOR. 
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In  the  first  portion  of  the  operate  cycle  diode  bridges 
in  Program  A  remain  non-conducting.  Hence  the  integrators 
will  operate  in  the  normal  fashion  until  x  =  L/2.  In 
program  B  all  bridges  are  conducting,  forcing  the  inputs 
to  all  integrators  to  zero.  These  integrators  will  remain 
holding  their  initial  conditions  until  x  =  L/2.  At  this 
instant,  the  comparator  will  switch  polarities  of  the  control 
voltages  to  the  bridges,  causing  those  in  program  A  to  conduct 
and  in  turn  forcing  the  input  to  each  of  these  integrators 
to  zero.  As  a  result  integrators  in  Program  A  will  hold 
their  value  from  x  =  L/2.  Switch  A  is  open,  hence  A  will 
have  no  effect  on  the  remainder  of  the  solution.  Diode 
bridges  on  integrators  in  program  B  will  now  be  non-conducting 
and  hence  program  B  will  now  be  in  the  operate  mode,  and  with 
switch  B  closed  will  generate  the  remaining  portion  of  the 
stored  function.  A  more  complete  description  of  the  operation 
and  performance  of  the  mode  control  of  the  integrators  used 
in  programs  A  and  B  is  given  in  Appendix  F. 

A  brief  description  of  the  operation  and  the  performance 
of  switches  A  and  B  will  now  be  given.  The  schematic  is 
given  in  Figure  lp-8.  Initially,  when  bridge  A  is  conducting, 
the  summing  junction  of  integrator  1  is  at  approximately 
0  volts.  If  the  protective  diodes  are  used  at  the  input  to 
the  diode  switch,  a  small  undesireable  step  occurs  at  the 
switching  instant  in  the  output  of  integrator  1.  If  how¬ 
ever,  the  diodes  are  left  out,  the  circuit  is  still  protected 
by  the  clamping  diodes  in  the  flip-flops  driving  the  bridge. 
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In  tnis  way,  switching  by  A  is  almost  ideal. 
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FIGURE  4-8.  SWITCHES  A  AND  B 


Diode  bridge  B  is  cut-off  in  the  first  portion  of  the 
operate  cycle.  To  prevent  the  input  side  of  the  diode  bridge 
from  rising  to  a  large  voltage  (to  the  input  voltage), 
protective  diodes  are  used  at  the  input  to  the  bridge. 

However,  at  the  instant  the  bridge  starts  to  conduct,  a 
small  undesireable  step  occurs  in  the  output  of  integrator  1 . 
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The  scheme  used  for  placing  a  diode  bridge  switch  in  an 
integrator  is  given  in  Figure  4-9. 


- — <  From  flip-flop 


o  Grid  of  Integrator 
(  0  volts) 


1  00K 


<From  flip-flop 

FIGURE  4-9.  DIODE  BRIDGE  SWITCH  IN  INTEGRATOR. 

When  the  bridge  is  cut  off,  the  diodes  have  a  small  charge 
stored  across  them.  At  the  Instant  the  bridge  switches,  the 
small  charge  stored  in  the  diodes  will  give  rise  to  a  small 
voltage  directly  to  the  grid  of  the  integrator.  This  will 
give  rise  to  a  step  of  approximately  0.3  volts,  the  sign  of 
which  is  dependent  on  the  input  voltage.  The  equivalent 
circuit  at  the  switching  instant  is  given  in  Figure  lj.-10. 


Small  charge  stored 
in  protective  diodes. 


oGrid 


FIGURE  4-10.  EQUIVALENT  CIRCUIT  FOR  DIODE  BRIDGE  SWITCH. 

An  improvement  by  a  factor  of  two  to  three  was  obtained 
using  the  circuit  shown  in  Figure  4*8* 
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4.3  FUNCTION  STORAGE  USING  NINE  POINTS 

The  program  for  the  storage  of  a  function  of  the 
independent  variable  x,  using  nine  points  will  now  be 
developed.  The  scheme  is  given  in  block  diagram  form 
in  Figure  Ll—  11. 


f (x, t ) 


O-Memorv 

Control 


X-Memory 
(nine  points) 

f (x, t ) 


FIGURE  1|-11  .  BLOCK  DIAGRAM  FOR  FUNCTION  STORAGE. 


4.3.1  THE  X-MEMORY 

The  X-Memory  is  analogous  to  the  sampler  previously 
described.  It  consists  of  a  series  of  nine  comparator 
controlled  memories  which  store  the  value  of  the  function 
at  equal  intervals. 

The  computer  scheme  for  the  input  to  the  X-Memory  is 
given  in  Figure  br-'\2 . 
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FIGURE  INPUT  TO  X-MEMORY. 


The  computer  scheme  for  the  Memory  storing  the  value  of 
the  function  at  x  =  L,  is  given  in  Figure  Integrator 

10  simply  generates  a  ramp,  the  value  of  which  is  proportional 
to  x.  Amplifier  44  will  track  the  input  voltage  until  x=  L, 
at  which  time  comparator  11  will  switch  the  mode  of  44  to  store. 
However  at  x  =  L,  the  computer  operate  cycle  is  completed  and 
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FIGURE  4-13.  X-MEMORY  STORING  f(L). 


the  computer  will  go  into  the  reset  mode.  This  would  mean 
that  comparator  11  would  immediately  switch  back  to  the  normal 
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state,  causing  the  X-Memory  to  again  track.  This  would 
not  provide  sufficient  time  for  the  0-Memory  to  read  in  the 
data.  For  this  reason  MM#3  is  used  to  provide  an  inhibit 
or  latching  voltage  to  insure  that  the  value  in  each  of  the 
X-Memories  is  held  for  10  milliseconds  after  the  operate  cycle 
is  completed. 

The  symbols  used  in  Table  lj.-2  listing  the  components  of 
the  X-Memory  are  explained  in  Figure  I4--II4.. 


(stored 
value ) 


C. A.  -  Comparator  Amplifier 

I.C.U.C.  -  Iterative  Control 
Unit  Comparator 

S.L.  Pot  -  Switching  Level 
Potentiometer 

S.A.  -  Storage  Amplifier 

D. B.  -  Diode  Bridge  for 

controlling  mode  of 
storage  amplifier 


C.T.  -  Capacitor  terminals 

for  storage  amplifier 

K  ,  Rj  ,  ...  denote  the  value 
of  the  function  at 
x  =  0,  L/8 ,  . . . 

I. A.  -  Input  amplifier 

(provides  function  to 
be  stored) 


FIGURE  [|-11+.  X-MEMORY  SYMBOLS 


39 


Value 

Stored 

S'. A. 

ICUC 

S.L. 

Pot 

S.L. 

Pot  S. 

D.B. 

C.T. 

C .  A . 

I. A. 

K0 

1+7 

1 

qi5 

.0915 

16 

B06 

BO  7 

15 

I.C. 

-  K-| 

48 

4 

<516 

.014.50 

17 

BO  8 

BO  9 

16 

P61 

+  K2 

49 

5 

Q2  5 

.1110 

18 

C06 

C07 

25> 

66 

-  k3 

32 

6 

Q26 

.1750 

19 

C08 

C09 

26 

F6i 

+  % 

33 

7 

Q35 

.2375 

30 

DO  6 

DO  7 

35 

66 

-  % 

34 

8 

Q36 

.2080 

21 

D08 

D09 

36 

p6l 

+  k6 

42 

9 

045 

.3642 

22 

10 

110 

45 

66 

-  K7 

43 

10 

046 

.4250 

23 

11 

111 

46 

^61 

+  K8 

44 

11 

Q4i 

.  4840 

29 

12 

112 

41 

66 

TABLE  14.-2.  X-MEMORY  FOR  NINE  POINT  FUNCTION  STORAGE. 

I4.  3. 2  O-MEMORY  CONTROL 

The  scheme  used  for  O-Memory  control  is  given  in  Figure 
I4-I9.  The  input  from  function  switch  I4I  to  amplifier  11 
provides  an  inhibit  pulse  when  a  scheme  for  obtaining  results 
on  an  x-y  plotter  is  used,  A  more  complete  explanation  for 
the  need  for  this  inhibit  pulse  will  follow  later. 

MM  wb  provides  an  inhibit  pulse  so  that  when  the  computer 
is  switched  from  the  “Initial  Condition  Read-In -position’’ 
to  "operate'*,  in  the  serial  solution  of  a  partial  differential 
equation,  at  least  one  iteration  for  the  solution  of  the  next 
interval  will  take  place.  This  will  insure  that  after  the 
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switch  from  I.C.  to  operate  is  made,  that  the  comparator  does 
not  switch  during  the  next  operate  cycle  causing  the  information 
stored  in  the  0-Memory  to  be  destroyed. 


diode  bridges  diode  bridges 

in  0-Memory  in  0-Memory 

FIGURE  4-1 5.  0-MEMORY  CONTROL. 


When  the  iteration  error  tends  to  zero,  that  is  when 
the  desired’ solution  has  been  obtained,  comparator  12  will 
switch  state  during  the  first  part  of  the  reset  cycle. 

This  will  cause  the  0-Memory  to  read  in  the  respective  values 
from  the  X- Memory  for  8  milliseconds.  The  X-Memory  has  been 
programmed  to  hold  its  value  for  10  milliseconds  in  the  reset 
cycle.  After  the  8  milliseconds  allowed  for  read  in,  an 
inhibit  pulse  is  obtained  from  MM  #2,  insuring  that  the 


0-Memory  will  hold  the  data  just  read  in  for  the  remainder 
of  the  reset  cycle  as  well  as  for  the  next  operate  cycle. 

Otherwise  comparator  12  would  remain  in  the  switched  state 

during  the  reset  cycle  as  well  as  the  next  operate  cycle* 

causing  the  data  read  in  from  the  X-Memory  to  be  lost. 

All  diode  bridges  in  the  0-Memory  are  switched  from 
amplifiers  £9  and  6£.  Potentiometer  Q^2  is  used  to  limit 
the  voltage  used  to  cut-off  the  diode  bridges.  Using  a 
smaller  voltage  to  reverse-bias  the  bridges  causes  a  marked 
improvement  in  reducing  drift  of  the  0-Memory  amplifiers  in 
the  store  mode. 

By  reading  in  data  for  8  milliseconds,  sufficient  time 
is  provided  to  use  larger  capacitors,  .01  ufd,  in  conjunction 
with  100K  resistors  in  the  0-Memories.  Error  during  read-in 
employing  these  values  is  approximately  0.U3  percent.  Using 
larger  capacitors  in  the  storage  amplifiers  of  the  0-Memory 
also  makes  the  real  time  plotting  scheme  more  feasible,  where 
the  value  in  the  0-Memory  must  be  held  for  at  least  seven 
seconds . 

14..3.3  0-MEM0RY 

The  0-Memory  consists  of  nine  memories  controlled  from 
Comparator  12  as  described  in  the  previous  section.  A  small 
bias  directly  to  the  grid  of  each  of  the  storing  amplifiers 
helps  to  overcome  drift  due  to  unbalance  in  the  diode  bridges. 
After  sufficient  time  has  been  allowed  for  computer  warm  up, 
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drift  in  the  0-Memory  can  be  limited  to  less  than  .02  volts 
for  10  seconds  using  this  technique.  Both  negative  and 
positive  compensation  has  been  provided.  The  schematic 
for  a  typical  0-Memory  is  given  in  Figure  L\.-'\6. 


(from  X-Memory) 


side  of  diode  bridge 
to  amplifier  5>9 
side  of  diode  bridge 
to  amplifier  65 


FIGURE  4-1 6.  TYPICAL  0-MEM0RY 


The  compensation  for  bridge  unbalance  is  given  in 
Figure  4-1  7- 
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FIGURE  4-1?.  COMPENSATION  FOR  DRIFT  IN  0-MEM0RY. 
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Table  4"  3  summarizes  the  amplifiers  used,  the  diode 
bridges,  capacitor  terminals,  as  well  as  the  required  drift 
compensation  for  the  0-Memory.  For  explanation  of  symbols, 
see  Figure  I4.-IU. 


Value 

Stored 

S  o  A . 

I. A. 

D.B, 

G.T. 

Drift  Compensation 

Comp . 

Pot . 

Pot . 
Setting 

+  Ko 

17 

47 

7 

00 

100 

A 

PI  7 

0  0865 

-  K1 

18 

48 

8 

01 

101 

B 

PI  8 

0 1  230 

+  k2 

19 

49 

q 

02 

102 

A 

PI  9 

.1390 

'  k3 

27 

32 

10 

0 

0  OJ 
Go 

A 

P27 

.1396 

+  K4 

28 

33 

1  1 

04 

104 

B 

P28 

.1200 

'  K5 

29 

34 

1  2 

105 

A 

P29 

.3380 

+ 

37 

42 

13 

06 

106 

B 

F37 

.1654 

-  K  7 

38 

43 

14 

~0T 

107 

A 

P38 

.0295 

+  kq 

39 

44 

15 

->■  0 
O  Gq 

GO 

A 

P39 

.0547 

TABLE  14.-3.  O-MEMORY  FOR  NINE  POINT  FUNCTION  STORAGE. 

4.3.4  NEWTON’S  ITERATIVE  FUNCTION  GENERATOR 

The  computer  scheme  used  for  programming  the  Newton’s 
Iterative  Function  Generator  is  given  in  Figures  i|-l8a, 

I4.-1 8b,  4-1 8c,  and  4-1 8d.  Integrator  91  generates  a  ramp 
which  is  proportional  to  x.  Comparator  performs  the 
switching  functions  described  by  Switches  A  and  B  while 
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comparator  #3  performs  the  function  of  holding  program  B  in 
initial  condition  for  the  first  portion  of  the  operate  cycle, 
then  switching  program  A  into  the  hold  mode  for  the  remaining 
half  of  the  operate  cycle.  The  arrangement  of  amplifiers 
90  and  81  is  necessary  since  the  removal  of  the  pot  set 
relay  from  integrator  5>U  means  that  the  initial  condition 
on  this  integrator  would  normally  be  -K0  -  f(0).  The 
portion  of  the  computer  designated  with  a  bar  above  all 
quantities,  for  example  e,  refers  to  program  B  of  Newton's 
Iterative  Function  Generator.  The  complete  computer  program 
for  the  function  generator  is  given  in  Figures  [j.-l8a,  I4.— 1 8b, 
ij.-l8c  and  i|-l8d. 


To  diode  switches  in  Automatic 
function  generator. 


FIGURE  14.- 1  8a. .  SWITCHING  FOR  THE  FUNCTION  GENERATOR 


FIGURE  J4.-1 8b . 
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FIGURE  JL|_— 1  8c .  ALGEBRAIC  CALCULATIONS  FOR 
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Figure  i8d,  computer  program  rai?  iterative  pumetionj  (3E/^tRAroR. 
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li.u  CONCULSIONS 

The  use  of  Newton’s  Iterative  Function  Generator  for 
extrapolating  a  stored  curve  provides  a  practical  means 
for  performing  function  storage  on  the  analog  computer., 

If  digital  storage  for  the  function  is  not  available, 
then  two  programs  of  five  points  represents  the  optimum 
choice  for  function  storage.,  The  reason  that  this  is  the 
optimum  choice  is  as  follows „  The  number  of  amplifiers 
required  for  two  programs  each  of  four  points  (total  number 
of  stored  points  is  seven)  is  twenty  amplifiers.  That 
for  two  programs  of  five  points,  making  a  total  of  nine 
stored  points  is  twenty-four.  It  is  evident  that  the  number 
of  amplifiers  required  per  point  for  the  two  programs  of  five 
points  is  less  than  that  for  the  two  programs  of  four  points. 
A  function  generator  consisting  of  two  programs,  each  of  five 
points  also  represents  an  optimum  choice  with  respect  to  the 
accuracy  of  the  function  generator.  In  the  calculations  for 
the  coefficients  a,  b,  c,  „.e,  gains  of  less  than  10  are  used 
on  all  amplifiers.  However  gains  of  twenty  are  required  on 
amplifiers  performing  these  algebraic  calculations  for  a 
seven  point  program  and  gains  of  thirty  five  are  required  on 
amplifiers  calculating  these  coefficients  in  a  program  for  8 
points  (see  Appendix  G) .  This  implies  that  if  the  stored 
function  is  a  ramp,  where  all  coefficients  following  b  are 
equal  to  zero,  the  errors  will  accumulate  rapidly  in  the 
algebraic  calculation  for  the  last  coefficients. 
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Pop  example,  in  the  algebraic  calculation  for  i,  the 
difference  between  eight  quantities,  each  with  large  gains 
(up  to  39  )  is  to  be  obtained.  This  will  tend  to  make  the 
required  potenti ometer  settings  extremely  critical. 

Using  field-effect  transistors  for  implementing  mode 
control,  Newton’s  Iterative  Method  may  be  used  for  parabolic 
curve  extrapolation  (that  is,  each  program  consists  of 
3  points)  using  any  number  of  stored  points.  The  iterative 
function  generator  program  however  will  only  require  fifteen 
amplifiers.  This  technique  makes  the  implementation  of 

function  storage  realizable  even  with  smaller  computers  such 
as  the  Pace  TR-I4.8. 

The  accuracy  of  the  extrapolation  of  the  stored  function 
performed  by  the  iterative  function  generator  is  dependent 
upon  the  irregularity  of  the  stored  function.  If  the  rate  of 
change  of  the  function  is  extreme,  then  more  points  will  be 
required  to  represent  the  function  more  accurately.  Using 
parabolic  extrapolation  with  mode  control  on  the  integrators 
performed  by  field-effect  transistors,  this  iterative  function 
generator  can  be  used  to  extrapolate  functions  with  only 
piecewise  continuous  derivatives.  The  error  of  the  extra¬ 
polated  function  at  the  instants  the  value  of  the  function 
is  stored  is  less  than  0.2  percent.  The  performance  of  the 

iterative  function  generator  for  nine  stored  points  is  given 
in  Figure  4-19. 
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TIME  (SECONDS) 


^.0  SOLUTION  OF  PARTIAL  DIFFERENTIAL  EQUATIONS  ON  THE 
ANALOG  COMPUTER 

In  this  chapter  the  solution  of  partial  differential 
equations  in  two  independent  variables  will  be  considered. 

There  are  three  general  techniques  that  can  be  applied  for 
the  analog  solution  of  these  equations. 

The  first  method  involves  treating  both  variables  in 
continuous  form.  Examples  include  rubber  sheet  and  soap 
film  analogies  (reference  13). 

A  second  approach  is  the  use  of  finite  difference  tech¬ 
niques  applied  to  both  independent  variables.  This  is  the 
same  technique  that  is  used  when  the  equation  is  being  solved 
on  a  digital  computer.  The  region  of  interest  is  divided 

into  a  mesh  and  a  system  of  equations  is  written  for  each 

2 

node.  The  total  number  of  nodes  is  N  if  the  range  of  each 
independent  variable  is  sub-divided  into  N  equal  intervals. 

The  application  of  this  scheme  is  not  feasible  on  the  analog 
computer  since  equipment  requirements  become  excessive. 

The  third  approach  to  the  solution  of  partial  differential 
equations  on  the  analog  computer  is  discreting  one  variable 
while  the  other  is  left  in  continuous  form.  This  method 
is  the  most  practical  of  the  three  described  for  utilization 
on  the  analog  computer. 

In  the  next  two  sections,  the  parallel  and  the  serial 
method  will  be  described.  In  each  of  the  methods,  finite 
difference  techniques  are  used  to  express  one  of  the 


' 
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independent  variables  in  discrete  form. 

5-1 .  PARALLEL  METHOD 

In  the  parallel  method,  the  range  of  the  x  variable 
is  divided  into  N  equal  intervals.  A  separate  computer  set¬ 
up  is  then  programmed  for  each  of  the  N-1  interpolation 
points.  The  solution  is  obtained  simultaneously  for  each  of 
the  interpolation  points.  The  advantage  of  this  method  is 
that  the  boundary  conditions  imposed  on  a  closed  region  may 
be  taken  into  account  automatically  without  requiring  any 
iteration  technique  to  find  a  remaining  boundary  condition 
on  the  solution.  The  unfavourable  aspect  of  using  this 
method  is  that  if  improved  accuracy  is  desired,  eauipment 
requirements  become  excessive.  The  general  scheme  involved 
in  the  parallel  solution  of  partial  differential  equations 
is  given  in  Figure  5-1  * 


FIGURE  5-1 .  PARALLEL  SOLUTION  OF  A  PARTIAL  DIFFERENTIAL 


EQUATION. 
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For  an  example,  consider  the  application  of  the  parallel 
method  to  the  solution  of  a  problem  in  heat  conduction. 

Problem 

A  rod  with  a  length  L,  with  isolated  lateral  boundaries 
has  the  end  x  =  0  maintained  at  T  =  T>|  when  t  >0,  while  the 
end  at  x  =  L  is  maintained  at  T  =  0,  The  inital  temperature 
distribution  along  the  rod  is  prescribed  as  T(x,u)  =  fix). 

The  equation  governing  the  distribution  of  heat  in  the  rod 
is  given  by; 


=  K  d  T  (x,  t )  (K  =  1  ) 
&t 


Solution 

Divide  the  length  L  into  N  equal  intervals.  The  length 
of  each  interval  is  then  L/N.  We  will  designate  x^  to  mean 
x  =  kL/N.  If  we  discrete  the  portion  dependent  upon  x; 


d  T(xn,t) 
dt 


T(xn+1,t)  -  2T(xn,t)  +  T(xn_>|,t) 

(ax)2 


If  the  equation  is  considered  at  station  x^  we  have; 

T(xp.t)  -  2T(xi ,t)  +  Ti  =  dT(xi ,t) 

(*x)2  dt 

The  computer  set-up  for  this  equation  at  station  x^  (which 
also  takes  into  account  the  boundary  condition  at  x  =  0)  is 
given  in  Figure  5>-2.  The  computer  may  be  programmed  in 
a  similar  fashion  for  stations  X£,  X3,  .  .  .  .xn_>j  . 
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FIGURE  5-2.  COMPUTER  SET-UP  FOR  STATION  X1  . 

5.2  THE  SERIAL  METHOD 

In  the  serial  method  of  computation,  the  portion  of 
the  partial  differential  equation  dependent  on  time  is 
discreted.  Then  in  the  calculation  for  the  solution  corres¬ 
ponding  to  t  +  ^t,  the  solution  for  time  t  must  be  available. 
That  is,  a  means  of  function  storage  must  be  available. 

The  difficulty  in  employing  this  method  is  that  if  the  partial 
differential  equation  is  being  solved  on  a  closed  region,  an 
iteration  scheme  must  be  devised  to  find  the  necessary  initial 
condition  at  x  =  0.  This  does  not  pose  a  very  serious  problem 
if  storage  is  available  on  the  analog  computer.  The  marked 
advantage  of  employing  this  method  is  that  the  equipment  is 
multiplexed.  The  same  computer  program  used  for  solving  the 
discreted  equations  at  time  t  is  also  used  for  the  next  time 
instant,  t  +  zxt.  Equipment  requirements  are  less  severe 
than  In  the  parallel  method.  By  decreasing  the  size  of  At 
more  solutions  will  be  obtained  in  a  given  region  of  interest, 
however  there  will  be  an  increase  in  any  inherent  instabilities. 


The  scheme  used  in  the  serial  solution  of  partial  differential 
equations  is  given  in  Figure  5-3*  The  arrows  indicate  the 
direction  of  integration. 
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FIGURE  5-3.  INTEGRATION  PROCEDURE  IN  THE  SERIAL  METHOD 


5-3.  SOLUTION  OF  A  HEAT  PROBLEM  USING  THE  SERIAL  METHOD 

The  serial  method  will  now  be  applied  to  a  problem 
in  heat  conduction  using  iterative  techniques  on  the  analog 
computer.  The  problem  that  will  be  considered  is  given 
below. 

A  rod  with  length  L,  with  isolated  lateral  boundaries 

has  ends  x  =  0  and  x  =  L  maintained  at  T  =  0  when  t  >  0. 

The  initial  temperature  distribution  along  the  rod  is 

prescribed  as  T(x,G)  =  5Q  sinfTX«  The  equation  governing  the 

L 

temperature  distribution  in  the  rod  is  given  by, 

~aT  ( x,  t )  =  K  ^2T(x,  t) 

}t 


K  =  L  =  1  (numerically ) 
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5*3.1  BLOCK  DIAGRAM  FOR  THE  ANALOG  COMPUTER 


If  we  now  discrete  the  portion  of  the  partial  differential 
equation  dependent  upon  time  we  have; 
d2T(x,tn)  +  d2T(x,tn_1)' 


dx ^  dx‘ 


where  tn  =  n^t 


=  T(x,tn)  -  T(x,  tn_;  )  - (5-1) 


a  t 


If  the  interval  At  is  made  small,  expression  (5-1 )  may  be 
approximated  by; 


d 

dx2 


T(x,tnj  -  T  ( x  ,  tn  - 1  )  . (?-2) 

At 


The  block  deagram  for  the  solution  of  equation  (5-2) 
using  iterative  techniques  on  the  analog  computer  is  given 
in  Figure  5-U*  The  required  function  storage  is  performed 
as  described  in  Chapter  l;.  The  computer  program  for  the 
X-Memory,  the  0-Memory,  their  respective  controls,  and  for 
Newton's  Iterative  Function  Generator  is  given  in  Chapter  !;. 


+T(0,tn) 


-T(x,tn_1 ) 


Program  for 
solving  equatiorf 
5-2. 


Mode  Control 
( I . C .  t o  Operate 


T(x, tn  ) 


Iteration 

Scheme 


1 


SPDT 

^  ....  .. 

Function 

Real  Time 

Relay 

Storage 
(see  Chap.!;) 

Plotting 

Scheme 

-<T(x,0)  from  DFG 


to  x-y  plotter 


FIGURE  5-lu  BLOCK  DIAGRAM  FOR  SOLUTION  OF  HEAT  PROBLEM. 
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5.3.2  COMPUTER  PROGRAM  FOR  SOLUTION  OF  EQUATION  (5-2) 

The  computer  program  for  the  solution  of  equation  (5-2) 
is  given  in  Figure  5-5. 


FIGURE  5-5.  COMPUTER  DIAGRAM  FOR  EQUATION  (5-2). 

Results  were  obtained  for  At  =  .0625.  The  potentiometer 
setting  used  for  this  particular  value  At  are  as  follows: 

Q00  -  .1000 
POO  -  .1600 
P01  -  .1600 

5.3.3  THE  ITERATION  SCHEME 

The  purpose  of  the  iteration  scheme  is  to  automatically 
derive  the  required  initial  condition,  T(0,tn).  The  computer 
program  is  given  in  Figure  5-6.  A  description  of  the  operation 
of  the  iteration  scheme  is  given  below. 

During  the  first  operate  cycle  where  the  temperature 
distribution  for  the  time  tn  is  being  sought,  an  initial 
guess  for  T'(0,tn)  is  provided  by  amplifier  21.  The  output 
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Figure  5*- 6.  puter  i  tef/? //on) 
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of  integrator  22  can  be  assumed  to  be  zero  because  of  the 
reset  operation  which  will  be  described  later.  The  input  to 
Amplifier  21  comes  from  the  armature  of  contacts  number  2  of 
relay  A,  which  is  driven  by  comparator  .  As  a  result,  the 
input  to  amplifier  21  is  an  initial  guess  adjusted  by  Ql4i  in 


the  case  where  T(0,  t)  is  being  sought,  or  (K>|  -  K0) 

Jr 


^n-1 


which  serves  as  a  good  approximation  for  T('0,tn).  Potentio¬ 
meter  Pipip  adjusts  the  initial  guess  so  that  it  is  slightly 
less  than  the  required  initial  condition,, 

Throughout  the  operate  cycle,  x-memory  23  will  track 
the  output  voltage  T(x,tnj.  At  x  =  L,  the  mode  of  x-memory 
23  will  change  and  the  value  of  T(L,t)  for  the  first  trial 
solution  will  be  stored.  This  value  will  be  designated  as 
T(L,tn)i#  Since  the  required  value  of  T(L,tn)  =  U,  T(L,tn)^ 
will  serve  as  an  error  for  the  first  iteration  and  will  be 
designated  as  „  O-memory  21;  will  now  track  held  by 

x-memory  23.  Read  in  can  take  place  for  the  first  ten 
milliseconds  of  the  reset  cycle,  because  the  mode  control 
for  x-memory  23  and  o-memory  21;  is  identical  to  the  X-Memory 
control  for  the  function  storage  (see  chapter  1;). 

During  the  reset  cycle,  diode  bridge  /f26  will  be  non¬ 
conducting,  hence  integrator  22  will  integrate  at  a  rate 
proportional  to  the  error  / ^  »  held  by  o-memory  23.  At  the 
end  of  the  reset  cycle,  diode  bridge  #26  will  again  conduct 
heavily,  causing  integrator  22  to  hold  its  last  value. 


Amplifier  21  will  now  provide  a  better  approximation  for 
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T(0,tn)  for  the  next  operate  cycle.  In  this  iterative 
fashion  the  correct  initial  condition  for  T(x,tn)  will  be 
obtained.  A  photograph  of  the  performance  of  the  iterating 
integrator  22  is  given  in  Figure  5-1 2b, 

Amplifier  u9  and  llj.  provide  mode  control  for  integrator 
22.  A  large  voltage  is  used  to  place  the  diode  bridge  into 
heavy  conduction  so  as  to  minimize  the  error  due  to  drift 
when  the  integrator  is  in  the  hold  state  (see  Appendix  F) . 
During  the  operate  mode  of  the  integrator  (when  bridge 
B  is  non-conducting),  a  smaller  voltage  is  used  to  place 
the  diode  bridge  into  the  non-conducting  state  in  order  to 
reduce  errors  due  to  unbalance  in  the  diode  bridge. 

After  k  iterations,  when  the  desired  solution  has  been 
obtained,  that  is,  T(L,tn)jc  =  0,  comparator  12  will  switch 
state.  This  will  cause  the  o-memory  in  the  function  storage 
to  read  in  the  values  of  T(x,tn).  The  comparator  will  remain 
in  this  state  for  8  milliseconds  at  which  time  MM  yl  will 
trigger  MM  #2.  This  will  place  an  inhibit  pulse  on  comparator 
12  insuring  that  this  comparator  switches  back  to  the  normal 
state  and  remains  in  this  state  for  the  remainder  of  the 
reset  cycle,  as  well  as  for  the  next  operate  cycle,  during 
which  an  iteration  is  started  for  obtaining  the  solution 
T(x,tn+1 ). 

At  the  instant  comparator  12  switches  state,  MM  ^6  is 
triggered,  producing  a  pulse  of  1  [}.  millisecond  duration. 
Amplifiers  2U  and  $6  place  diode  bridge  C  into  the  conducting 
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state  for  the  duration  of  the  pulse.  This  will  cause 
integrator  22  to  reset  to  zero  since  the  output  voltage 
from  o-memory  2l\.  is  zero  when  the  desired  solution  has 
been  obtained. 

The  output  voltage  from  amplifier  21  is  now  equal  to 
( Ki  -  K0)^n  and  the  output  from  Newton's  Iterative  Function 
Generator  is  T(x,tn).  The  computer  will  now  iterate  for  the 
solution  T(x,tn+^)  in  a  manner  analogous  to  that  described 
for  T(x, tn) . 


5.3.i|  MODE  CONTROL  (I.C.  TO  OPERATE) 

Mode  control  is  required  in  the  solution  of  the  heat 
problem  to  change  the  input  (corresponding  to  T(x,tn_^)  in 
Figure  5“5)  to  the  first  integrator  of  the  system  solving 
Equation  (5-2).  This  input  is  changed  from  the  diode 
function  generator  output  which  generates  the  temperature 
distribution  at  t  =  0,  to  the  output  from  Newton's  Iterative 
Function  Generator.  That  is,  when  the  solution  for  T(x,Atj 
has  been  obtained,  the  mode  of  the  system  must  be  changed 
so  that  the  remaining  solutions  T(x,kAt)  may  be  obtained. 

The  computer  program  for  this  mode  control  is  given  in 
Figure  5“7>  followed  by  a  description  of  its  operation. 
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FIGURE  5-7.  COMPUTER  PROGRAM  FOR  MODE  CONTROL. 


The  operation  of  the  mode  control  can  be  best  explained 
by  considering  the  input  voltage  to  the  comparator  Nq .  This 
is  depicted  graphically  in  Figure  5-8. 

If  function  switch  11  is  closed,  the  input  to  comparator 
N-j  will  always  be  negative.  This  comparator  will  remain  in 
what  will  be  referred  to  as  the  switched  state.  The  input  to 
integrator  00  (see  Figure  5-5)  will  be  -T(x,0),  the  initial 
temperature  distribution  in  the  rod.  When  the  iterations 
for  the  solution  of  T(x,Atj  are  completed,  the  iteration 
scheme  will  reset  and  the  resulting  function  T(x,At)  will  be 
read  into  the  0-Memory  of  Newton's  Iterative  Function 
Generator.  However,  since  the  input  to  integrator  00  will 
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FIGURE  5-8.  INPUT  VOLTAGE  TO  COMPARATOR  . 

still  remain  T(x,0),  this  same  solution  will  be  continually 
solved  until  the  state  of  relay  A  is  changed. 

To  obtain  the  solution  T(x,2At),  the  state  of  comparator 
N>j  must  be  changed  so  that  when  T(x,etj  is  obtained,  the  state 
of  relay  A  will  be  such  that  -T(x,At)  will  be  feed  into 
integrator  uu  for  the  next  operate  cycle.  However,  this 
transition  must  occur  only  during  the  reset  cycle  if  the  data 
stored  in  the  0-Memory  is  not  to  be  destroyed.  The  logic 
required  for  accomplishing  this  operation  is  performed  by 
amplifier  30  and  comparator  N-j  .  Function  switch  11  may  be 
opened  while  the  computer  is  either  in  the  operate  or  the  reset 
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mode.  These  two  cases  are  now  investigated  using  Figure 

5-8. 

In  the  first  case,  consider  FS  1 1  being  opened  when  the 
computer  is  in  the  operate  mode.  The  voltage  to  comparator 
N'l  will  rise  to  -0.5  volts.  The  state  of  the  comparator  N^ 
will  not  change.  At  the  instant  the  computer  switches  into 
the  reset  mode,  the  input  voltage  to  comparator  will  rise 
to  1.5  volts.  This  will  cause  to  switch  state,  driving 
relay  A  into  its  normal  state.  The  switching  speed  of  the 
relay  used  is  in  the  order  of  five  milliseconds,  therefore 
there  will  be  sufficient  time  in  the  reset  cycle  for  the 
relay  to  change  its  state.  When  the  computer  returns  to  the 
operate  mode  ( chopper  drive  voltage  becomes  -I4.  volts)  the 
input  voltage  to  comparator  becomes  -0.5  volts.  However, 
because  of  the  one  volt  dead-zone,  comparator  will  remain 
in  its  normal  state. 

In  the  second  case,  FS  1 1  is  opened  when  the  computer 
is  in  the  reset  mode.  No  complications  arise  in  this  case. 
Comparator  will  switch  to  the  normal  state  and  remain  in 
this  state  as  described  in  the  previous  paragraph. 

When  FS  1 1  is  opened,  no  information  is  available  as  to 
what  state  the  iteration  scheme  is  in.  However,  within 
0.1  seconds  the  first  solution  T(x,At)  is  obtained  and  is 
stored  in  the  0-Memory.  At  the  instant  comparator  switches 
state,  MM  y6  is  triggered,  resetting  integrator  22.  An  in¬ 
hibit  pulse  is  also  placed  on  comparator  12  by  MM  insuring 
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that  least  one  more  iteration  is  performed  before  this 
comparator  will  switch  state.  The  input  to  integrator  00 
is  now  equal  to  T(x,^t),  which  is  obtained  from  the  Newton's 
Iterative  Function  Generator.  The  computer  will  begin  itera¬ 
tions  for  T(x,2/^t).  The  remaining  solutions  T(x,kAt)  are 
obtained  in  a  manner  described  in  Section  5.3.3. 

5.3.5  REAL  TIME  PLOTTING  SCHEME 

This  scheme  was  devised  in  order  to  obtain  results 
T(x,k4t)  on  an  x-y  plotter.  The  essence  of  the  scheme  is  as 
follows : 

When  a  solution  for  a  new  increment  of  time  is  obtained, 
this  solution  is  first  read  into  the  0-Memory  of  Newton's 
Iterative  Function  Generator.  An  inhibit  pulse  is  then  placed 
on  Integrator  22  so  that  no  convergence  to  the  solution  for 
the  next  increment  of  time  is  obtained.  A  sampled  data 
plotting  scheme  similar  to  that  described  in  reference  11  is 
then  used  to  read  out  the  function  stored  in  the  0-Memory. 

Five  seconds  are  allowed  for  the  x-y  plotter  to  record  the 
solution  and  two  seconds  for  the  reset  of  the  plotter.  The 
inhibit  pulse  on  integrator  22  is  then  removed  while  the  com¬ 
puter  is  in  the  reset  mode  and  iterations  for  the  solution 
for  the  next  interval  of  time  will  begin.  The  entire  solution 
T(x,kAt)  may  be  obtains  d  on  an  x-y  plotter  by  a  repetition  of 
the  above  process.  The  computer  program  for  the  real  time 
plotting  scheme  is  given  in  Figure  5-9.  The  prefix  E  will 
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FIGURE  ^-9.  HEAL  TIME  PLOTTING  SCHEME. 
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imply  that  the  equipment  designated  in  this  way  is  part  of 
the  Pace  TR-1U  computer.  A  more  detailed  description  of 
the  operation  of  the  plotting  scheme  is  given  below. 

When  the  iteration  for  a  particular  solution  is 
completed,  comparator  12  will  switch  state,  triggering  MM  #1. 
Bight  milliseconds  later  MM  jf2  will  be  triggered.  This  will 
cause  comparator  N2  to  change  state,  driving  relay  B  to  the 
switched  state.  Integrator  E7  will  now  produce  a  ramp  which 
will  reach  ten  volts  in  five  seconds.  This  output  serves 
as  a  drive  for  the  x-axis  of  the  plotter  and  as  an  input  for 
the  sampled-data  plotting  scheme  comparator  08.  As  the 
computer  produces  the  solution  in  repetitive  operation, 
x-memory  E9  will  sample  the  function  T(x,k^t)  at  successive 
intervals  of  x.  The  o-memory  will  hold  or  store  these  values 
and  serve  as  the  input  to  the  y-axis  of  the  plotter.  In  five 
seconds,  1 ,430  samples  of  the  function  will  be  taken.  Thus, 
the  output  of  o-memory  El  0  will  approach  a  continuous  curve. 

At  the  instant  comparator  N 2  switches  state,  amplifier 
E6  will  place  an  inhibit  pulse  on  integrator  22  causing  it 
to  reset  to  zero  and  remain  in  this  state  until  the  plot  on 
the  x-y  plotter  is  completed.  At  the  instant  Integrator  E7 
reaches  10  volts,  comparator  07  will  change  state  causing 
comparator  Np  to  return  to  the  normal  state.  This  will  cause 
integrator  E7  to  reset  to  zero  volts  with  a  time  constant 
equal  to  .1  seconds. 

The  input  voltage  from  amplifier  Elf.  will  cause  comparator 
H,  to  provide  an  inhibit  pulse  on  Integrator  22  for  an 
additional  two  seconds.  The  operation  of  comparator 
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Nj  will  be  explained  by  reference  to  Figure  5_/tU,  which 
gives  the  input  voltage  to  this  comparator. 


FIGURE  5-10.  INPUT  VOLTAGE  TO  COMPARATOR  . 

The  reference  voltage  to  comparator  N^  at  the  instant 
comparator  N^  switches  back  to  the  normal  state  is  -7.5  volts. 
This  reference  voltage  will  decay  to  +.5  volts  in  2  seconds. 
This  repons e  is  obtained  from  the  program  of  a  damped  second 
order  system.  Superimposed  on  the  chopper  d^ive  is  a  pulse 
of  5  millisecond  width  which  is  derived  from  MM  #5.  The 
purpose  of  this  pulse  is  to  insure  that  when  the  Inhibit 
pulse  is  removed,  it  is  in  the  first  5  milliseconds  of  the 
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reset  cycle.  Otherwise  if  switching  occurs  in  the  latter 
part  of  the  reset  cvcle ,  a  risk  of  destroying  the  stored 
function  exists.  When  comparator  switches  back  to  its 
normal  state  it  will  remain  in  this  state  because  of  the 
1  volt  dead  zone. 

When  the  inhibit  pulse  from  comparator  is  removed, 
the  iteration  scheme  will  proceed  to  find  the  solution  for 
the  next  interval  of  time  in  the  normal  fashion.  When  the 
next  solution  is  obtained, the  sequence  of  operations  for 
plotting  this  solution  on  an  x-y  plotter  is  identical  to 
that  for  the  previous  solution. 

Amplifier  E5>  and  comparator  now  perform  the  mode 
control  described  in  Section  5.3.U.  With  S>|  closed, 
comparator  remains  in  the  switched  state.  The  computer 
will  continue  obtaining  the  solution  T(x,At).  By  use  of 
PS  30  both  T(x,0)  and  T(x,At)  may  be  plotted.  With 
open,  no  change  in  the  state  of  comparator  N>j  will  occur 
until  the  inhibit  pulse  from  comparator  N3  is  removed.  A 
better  explanation  for  the  requirement  that  Nj  must  switch 
in  the  first  5  milliseconds  of  the  reset  cycle  may  be  given. 

The  switching  time  for  relay  A  driven  by  comparator  is  in 
the  order  of  5  milliseconds.  Therefore  a  safety  factor  of 
5  milliseconds  exists  for  relay  A  to  completely  switch  in  the 
reset  cycle.  If  switching  is  not  completed  in  the  reset  cycle 
tte  solution  T(x,At)  could  be  destroyed,  rendering  the 
remainder  of  the  solutions  meaningless. 
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5.3.6  INITIAL  CONDITION  T(x,0) 


The  initial  temperature  distribution  in  the  rod  is 
generated  using  a  diode  function  generator.  The  curve  is 
approximated  using  16  segments „  The  computer  program  is 
given  in  Figure  5-1 1 . 

+100 


o 


00 

FIGURE  5-11.  GENERATION  OF  INITIAL  CONDITION  T(x,0). 


5.4  THEORETICAL  RESULTS 

The  separation  of  variables  method  will  be  applied  to 
the  following  partial  differential  equation. 

DT(x,t)  =  DT(xtt)  0  £  x  ^L  =  1 

with  T(x,0)  =  50  sin  x/L 
T(0,t)  =  0 
Tf L, t )  =  0 

Assume  a  solution  of  the  form;  T(x,t)  =  T(t)XCx).  Substituting 
this  in  the  partial  differential  equation  we  have; 

X"T  =  XT  ' 
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T7e  now  have  two  ordinary  differential  equations.  The  first 
one  i s ; 

T  +  72T  =  0 

-A2t 

therefore  T(t)  =  c  e  A 

1 

The  solution  for  the  space  dependent  part  is 

X"  +  X2X  =  0 

therefore  X  =  cpsin2x  +  c^cos^x 

If  the  initial  conditions  are  now  applied: 

X(0)  =  0  therefore  c^  =  0 

X(L)  =  0  =  C2sin^L 

therefore  -a L  =  nrr 

L 

The  complete  solution  is; 

T(x,t)  = 

T(x,0)  = 


cn  gin  nrrx  e 


2  2 

-n  rT  t 


IT 


cR  sin  nfrx 


But  this  is  just  the  Fourier  Series  for  T(x,0)  where  the 
coefficients  cn  are  given  by 

c. 


n 


=  2  \  T(x,0)  sin  nrrx  dx 

L  JO  L 


Substituting  L  =  1,  and  T(x,0)  =  5>0  sinrrx, 

/I 


n 


=  2 


0 


5>0  sinrrx  sin  nrrx  dx 


The  integral  is  non-zero  iff  n  =  1 ,  at  which  time 

=  1 00  \  sin 2rfx  dx  =  100^  (1  -  cos  2rfx)  dx  =  £0 
JO  Jo  2 


'1 


The  theoretical  solution  is, 
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2 

T(x,t)  =  50  e  ^  ^  sinrrx  . (5-3) 


5*5  COMPARISON  OP  THEORETICAL  AND  COMPUTER  RESULTS 


Photographs  of  results  obtained  on  the  analog  computer 
are  given  in  Figure  5-1 2c.  The  smallest  value  of  At  that 
could  be  chosen  was  .0625.  The  reason  for  this  restriction 
will  be  given  after  the  theoretical  solution  is  obtained 
for  the  difference  equation.  Consider  the  following  discreted 
equation  which  gives  the  temperature  distribution  at  time  t. 

. (5-Ji) 


d2T(x,t)  =  T(xft;  -  T ( x , t — At ) 

2 


dx 


a  t 


Let  l/^t  =  a 

The  theoretical  solution  obtained  for  this  problem  showed 
that  the  temperature  distribution  for  any  time  t  was  of  the 
form,  B  sinux.  Therefore  we  assume  that, 

T(x,t-At)  =  B  sinrTx 
T( 0 , t )  =  A 

Substituting  the  above  quantities  into  equation  C5-lx)  we 


have , 

d^T(x,  t )  ~  a  T(x,t)  =  -a  B  sinrTx  . . 

dx2 

If  we  now  take  the  Lanlace  transform  of  equation  (5-5); 
S2t(s)  -  sT(0+ )  -  T(0+)  -  aT(s)  =  -aBrr 


(5-5) 


S2  +  FT2 


( s 2  -  a)T(s)  =  A  -  aBrr 


s2  +  rr2 


T(s)  = 


aB  it 


(s*  -  a)  (s' 


a 


)  (  s2  +  it 2) 
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Using  Newton’s  partial  fraction  expansion  method  we  obtain; 


T  ( s )  =  aB  rr  -  aBrf  1 

( IT  2  +  a )  C  s  -  +  rr  2 )  2  ifa  ( a  +  rr  2 )  (  s  -  'Ta ) 

+  aB  tt _  _ 1  +  A  _ 1 _ 

2iTa(a  +  rT2)  (s  +  J~a)  2  JIT  (s  -  fa) 

-  _a_  1  . (5-6) 

(  s  +  \fa ) 


The  corresponding  time  solution  for  equation  (5-6)  is; 


T(x,t)  = 
+ 

where 


■  ^  B  sinrrx  -  C2  B  e 


c^  A  e 


■'Ja  ; 


C3  A  e 


+  (Ta  x 
-fa  x 


C2  B  e~^  x 


(5-7) 


a 


Crr^  +  a) 
a  iT 


'3 


(2  JliHa  +  1T2) 

1 


2\ja 


In  order  that  the  boundary  condtion  at  x  =  1  be  satisfied, 
(that  is,  T ( 1 , t )  =  0  ); 

C2B  =  c^A 


A  =  c2B 

c3 

Equation  ($-7)  will  be  used  to  calculate  the  theoretical 
solution  to  the  discrete  equation  .  The  actual  solution 

to  the  problem  will  be  obtained  using  equation  (5-3).  Table 
5-1  gives  a  comparison  between  the  computer  solution,  the 
theoretical  solution  to  the  discrete  equation  and  the  actual 
solution  for  a  time  increment  of  .0625  seconds.  In  Table 
5-1,  M  will  signify  the  maximum  value  of  the  sinusoidal 
solution  for  each  interval  of  time. 
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Time 
(sec . ) 

Computer  Results 

Theoretical  Soln 

(using  eqn.5-7) 

M 

Theoretical  Soln 
(using  eqn.  5-3) 
M 

0 

50 

50 

50 

.062$ 

31 

30.9 

26.9 

•  12$0 

19 

19.1 

1ii.$ 

.187$ 

12 

1 1  .8 

7.8 

.  2$00 

7 

7.3 

k.2 

TABLE  5-1.  COMPUTER  AND  THEORETICAL  RESULTS. 

In  order  that  the  solution  of  the  discreted  equation 
(5-4)  give  a  closer  approximation  to  the  actual  solution, 
a  smaller  increment  in  time  must  be  chosen.  The  computer 
solution  of  the  discreted  equation  compares  very  favorably 
with  the  theoretical  solution  obtained  using  equation  (5-7). 
Using  the  real  time  plotting  scheme  the  error  can  be  shown  to 
be  less  than  .5  percent.  A  comparison  between  the  solution 
of  equation  (5-4)  an^  the  actual  solution  given  by  equation 
(5-3)  is  given  in  Table  5-2  for  At  =  .01  seconds. 

Prom  'fable  5-2,  the  error  using  the  discrete  equation 
(5-4;  is  2.4  percent  as  compared  to  1 5  percent  for  approxi¬ 
mately  the  same  time  instant  using  At  =  .0625  (see  Table 
5-1).  If  the  increment  in  time  is  made  smaller,  even  better 
results  may  be  obtained. 

An  explanation  of  the  reason  why  the  increment  in  time 
could  not  be  made  less  than  .06  seconds  will  now  be  given. 
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Time 

( seconds ) 

Theoretical  Solution 
(using  equation  5-7) 

M 

Theoretical  Solution 
(using  equation  5-3) 

M 

0 

50 

50 

.01 

4?. 5 

47.1 

.02 

111  .k 

41 .1 

.03 

37.7 

37.2 

.04 

34.3 

33.6 

.05 

31 .2 

30.5 

.  06 

28.4 

27.6 

TABLE  5-2.  THEORETICAL  RESULTS  FOR  At  =  .01  . 

A  careful  examination  of  equation  (5-7)  will  show  that  if 
the  value  of  A  is  corrupted  by  noise  of  a  magnitude  S ,  then 
the  solution  will  deviate  from  the  desired  solution  by  £e^  x. 
At  the  final  value  of  x,  namely  x  =  1,  the  solution  will 
diverge  by  £ea.  For  the  case  of  At  =  .0625,  this  deviation 
is  equal  to  6kS •  An  examination  of  the  initial  condition 
on  integrator  00  showed  a  maximum  possible  error  of  .2  volts. 
This  offers  sufficient  explanation  why  the  solutions  tend 
to  ."’flutter4  towards  the  end  of  the  operate  cycle  when  the 
computer  is  in  repetitive  operation.  By  making  At  even 
smaller,  this  error  is  excentuated  even  further. 

The  source  of  this  error  was  traced  to  the  choppers 
stabilizing  the  amplifiers  in  the  computer.  Mode  control 
for  the  integrators  was  then  built  up  using  field-effect 
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transistors  switched  from  a  General  Radio  pulser.  This 
pulser  was  triggered  externally  from  the  9  it  cycle/sec. 
voltage  driving  the  stabilization  choppers.  Results  obtained 
using  this  external  mode  control  showed  no  perceptible 
"flutter" . 

It  may  be  concluded  that  the  application  of  iterative 
techniques  to  the  solution  of  partial  differential  equations 
on  the  analog  computer  represents  a  powerful  method. 

Examining  Table  %-Z  shows  that  the  error  between  the 
theoretical  and  the  computer  solution  to  the  discrete 
equation  (5>-lf)  is  almost  insignificant.  With  external 
mode  control  implemented  with  field-effect  transistors,  much 
smaller  values  for  the  increment  in  time  (At)  may  be  used, 
giving  results  well  within  one  percent  error.  Using 
transistorized  operational  amplifiers  (Nexus  FSL-12) 
with  mode  control  performed  by  field-effect  transistors, 
iterative  computation  may  easily  be  performed  at  10  kilocycles 
per  second  (reference  1i|).  Since  these  operational  amplifiers 
use  no  choppers  for  stabilization,  the  error  previously 
described  will  no  longer  exist.  This  would  then  present  a 
practical  technique  for  the  solution  of  partial  differential 
eauations  using  iterative  techniques  at  ultra-high  speed  on 
the  analog  computer. 
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FIGURE  5-1 2a.  ITERATIVE  SOLUTION  OF  BOUNDARY  VALUE  IN  O.D.E. 


FIGURE  5-1 2b.  PERFORMANCE  OF  ITERATING  SCHEME. 
(INPUT  TO  COMPARATOR  12) 


FIGURE  5-1 2c .  SOLUTION  OF  HEAT  PROBLEM  (^t  =  .0625). 
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( 1  ip )  ASH,  J.A:  The  use  of  field-effect  transistors  for 

implementing  mode  control  for  analog  computers 
was  demonstrated  in  the  Electrical  Engineering 
Department  of  the  University  of  Alberta. 


APPENDIX  A:  ITERATIVE  CONTROL  UNIT 


The  Iterative  Control  Unit  consists  of  the  specialized 
equipment  required  for  iterative  computation  on  the  analog 
computer.  The  unit  includes  twelve  electronic  comparators, 
three  comparators  with  variable  dead  zone,  thirty- six  diode 
bridges,  six  monostable  multivibrators,  and  one  double 
pole  double  throw  relay.  This  Iterative  Control  Unit  also 
contains  its  own  plus  and  minus  twelve  volt  power  supply 
(Dressen-Barnes  Electronic  Corporation  Power  Supply, 
model  20-1 2S).  A  front  and  side  view  of  the  Iterative  Control 
Unit  giving  the  equipment  layout  is  given  in  Figure  A-1 . 


FIGURE  A-1.  ITERATIVE  CONTROL  UNIT. 
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The  front  panel  of  this  unit  has  been  drilled  to  accept 
the  shielded  leads  used  in  patching  on  the  Face  Computer. 

This  faciliates  direct  patching  between  the  computer  patch 
panel  and  the  Iterative  Control  Unit.  Details  of  the  method 
by  which  contact  is  made  from  the  patch  panel  on  the  computer 
to  the  Iterative  Control  Unit  is  given  in  Figure  A-2. 


outer  insulation 
shield 


gold-plated 

contac  t 


front  panel 
Iterative 
Control  Unit  (drilled 
to  accept  patch  cords 
for  computer) 


printed 

circuit 

"board 


berilium 
copper  wire 


Berilium 
copper  wire 
is  soldered 
inside  turret 
terminal 


FIGURE  A-2.  CONTACT  ARRANGEMENT. 


A  more  detailed  description  of  the  individual 
components  of  the  Iterative  Control  Unit  is  given  in  the 
Appendices  that  follow. 


. 


APPENDIX  B.  DIODE  BRIDGE  SWITCH 


The  major  component  in  implementing  storage  on  an 
operational  amplifier  is  the  diode  bridge  switch.  The  need 
for  an  ideal  electronic  switch  poses  a  serious  problem  in 
accurate  iterative  computation.  In  the  diode  bridges 
incorporated  into  the  Iterative  Control  Unit,  all  resistors 
were  matched  within  .1  percent  while  the  forward  drop  across 
the  diodes  in  conduction  was  matched  within  .5  millivolts. 

A  schematic  for  a  diode  bridge  switch  is  given  in  Figure  B-1  . 


All  diodes 
HS  1020. 


FIGURE  B-1.  DIODE  BRIDGE  SWITCH  SCHEMATIC. 


The  scheme  used  for  measuring  the  forward  drop  across 
the  diodes  in  conduction  is  given  in  Figure  B-2. 

When  used  in  a  track-hold  circuit,  the  error  due  to  the 
unbalance  in  the  diode  bridge  is  less  than  0.0£  volts. 
However  in  the  case  where  a  diode  bridge  is  used  in  the 
summing  junction  of  an  integrator  (Switch  A  and  B  in  Newton* 
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+1 2  volts 


Differential  AC-DC 


FIGURE  B-2 .  MEASURING  OF  FORWARD  RESISTANCE  OF  DIODE. 


Iterative  Function  Generator;  an  unbalance  of  .u5>  volts 
integrated  for  ten  seconds  gives  an  error  of  .5  volts. 

To  reduce  the  unbalance  in  these  particular  diode  bridges 
even  further,  a  modification  to  the  diode  bridge  was  made  as 
indicated  in  Figure  B-3. 


The  bridge  can  be  nulled  almost  exactly  giving  an 
integrated  error  for  ten  seconds  of  less  than  0.03  volts. 
Better  accuracy  could  be  obtained  if  all  diode  bridges  were 
placed  in  an  oven  thereby  reducing  drift  due  to  changes  in 
temperature. 


APPENDIX  C.  ELECTRONIC  COMPARATOR 


In  order  to  achieve  faster  switching  for  the  diode 
bridges,  an  electronic  comparator  was  built  up.  A  block 
diagram  for  the  comparator  is  given  in  Figure  C-1 . 


input 


+11 .35N 
-U.6  S 
-0.6  N 
+11. 35s 


-ii .35n 
+u  .6  S 

+u  o  6  N 
-11 .35s 


FIGURE  C-1.  COMPARATOR  BLOCK  DIAGRAM. 

The  symbols  N  and  S  refer  to  the  normal  and  switched 

states  respectively  of  the  comparator.  The  performance  of 

each  of  the  components  of  the  comparator  are  as  follows: 

Schmitt  Trigger:  Rise  time  -  .  2/^seconds 

Fall  time  -  .2/+ seconds 

Trigger  point  1 .54  volts  (see 

Table  C-1  for  the 
trigger  point  of  ea 
comparator) 

An  emitter  follower  input  stage  with  a  clamping  diode  on  the 
input  allows  the  input  voltage  to  the  comparator  to  swing 
from  -200  to  +200  volts. 

Direct  coupling  is  used  from  the  Schmitt  trigger  to  the 
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bistable  multivibrators.  This  is  done  to  insure  the  correct 
state  of  each  of  the  multivibrators.  Clamping  diodes  are  used 
on  the  collectors  of  the  transistors  to  prevent  the  excursion 
into  the  reverse  biased  region  from  exceeding  approximately 
i  volt.  This  negative  excursion  into  the  reverse  biased 
region  is  required  to  insure  that  the  diode  bridges  driven 
by  the  bistable  multivibrators  are  reverse  biased  or  non¬ 
conducting.  Otherwise  if  the  collector  voltage  of  the  con¬ 
ducting  transistor  approached  zero  volts  (or  a  slightly 
positive  voltage)  this  x«rould  not  insure  that  this  voltage 
applied  to  the  diode  bridges  would  stop  them  from  completely 
conducting.  The  performance  of  the  bistable  multivibrators 
is  as  follows: 

Rise  time  =  0 .ipuseconds 

Pall  time  =  0 . 7t^seconds 

A  test  to  check  the  performance  or  the  balance  of  the  voltage 
outputs  from  the  n-p-n  and  p-n-p  bistable  multivibrators  is 
given  in  Figure  C-2. 


FIGURE  C-2.  PERFORMANCE  TEST  FOR  COMPARATOR. 


-  86 


As  shown  in  Figure  C-2,  the  unbalance  gives  rise  to  a 
spike  of  six  to  seven  volts  with  a  duration  of  less  than 
0. 2^ seconds .  This  switching  effect  causes  an  error  less 
than  0.05  volts  in  a  memory  using  a  .001  ufd.  capacitor. 

The  switching  level  for  each  comparator  is  given  in  Table  C-1 . 
The  wiring  diagram  for  the  comparator  is  given  in  Figure  C-3. 


Comparator 

Switching 

Level 

Comparator 

Switching 

Level 

1 

0.85 

7 

1 .25 

2 

1 .58 

8 

1.14.5 

3 

1 .60 

9 

1  .08 

u 

1 .75 

10 

1 .25 

£ 

1 .40 

1 1 

1.30 

6 

1 .25 

1  2 

1 .214. 

TABLE  C-1  .  SWITCHING-  LEVELS  OF  COMPARATORS 
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All  diodes 


+  1  2N 
-0.63 

+  1  23 
-0.6N 


-123 

+0 . 6N 


-1  2N 

+0.6S 


2N1 309 
2N1 308 
T-TS1020 


FIGURE  C-3.  WIRING  DIAGRAM  FOR  COMPARATOR. 


APPENDIX  D.  NEXUS  COMPARATOR 


The  need  for  a  fast  comparator  with  a  variable  dead 
zone  arises  in  the  application  of  syncronization  techniques 
as  well  as  in  the  real  time  read  out  scheme  in  iterative 
analog  computation.  The  basis  for  this  comparator  is  a  fast 
transistorized  operational  amplifier,  the  Nexus  PSL  -  6  or 
the  FSL-12.  Provision  was  made  for  this  operational  amplifier 
to  be  used  either  as  a  unity-gain  inverter  or  as  a  variable 
dead-zone  comparator  with  limiters  on  the  output. 

When  used  as  an  inverter,  the  circuit  will  work  well 
up  to  1  (JO  kilocycles.  When  used  as  a  comparator,  the  rise 
and  fall  times  are  in  the  vicinity  of  four  to  five  micro¬ 
seconds  with  the  limiters  set  to  limit  the  output  voltage  to 
plus  and  minus  five  volts.  With  a  dead  zone  of  one  volt, 
this  circuit  will  operate  in  a  realiable  fashion  to  about 
forty  kilocycles.  The  wiring  diagram  for  the  switch  in  the 
aInverterir  and  in  the  ‘‘'Comparator1'  position  is  given  in 
Figure  D-1  and  Figure  D-2  respectively. 


1  OK 

e^n  o - v'Vw— 


eref  ° 


1  UK 
vwv- 


1  OK 
www — 


— c 

— 

+ 

/^FSL- 

j>—  - 

//  (FSL- 

out 


All  resistors  1  watt 


FIGURE  D-1  . 


NEXUS  OPERATIONAL  AMPLIFIER  USED  AS  AN  INVERTER. 
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FIGURE  D-2.  WIRING  DIAGRAM  FOR  VARIABLE  DEAD-ZONE  COMPARATOR. 


APPENDIX  E.  MONOSTABLE  MULTIVIBRATOR 


The  purpose  of  the  monostable  multivibrator  is  to  provide 
an  inhibit  pulse  of  fixed  duration  in  iterative  analog  compu¬ 
tation.  Figure  E-1  gives  the  wiring  diagram  for  a  monostable 
multivibrator  which  will  give  a  negative  inhibit  pulse. 

Figure  E-2  gives  the  wiring  diagram  for  one  that  will  give 
a  positive  pulse.  The  width  or  duration  of  the  inhibit  pulse 
in  each  case  is  given  approximately  by; 

3 

T  =  82K  x  C  x  ln2  =  10  seconds. 


FIGURE  E-1.  MONOSTABLE  MULTIVIBRATOR  FOR  NEGATIVE 

INHIBIT  PULSE. 


INHIBIT  PULSE. 
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Table  E-1  gives  a  tabulation  of  the  various  monostable 
multivibrators  used  in  the  Iterative  Control  Unit.  Duration 
of  the  inhibit  pulse  and  the  voltage  output  in  the  normal 
state  for  each  is  given. 


Monostable 

Multivibrator 

Sign  of 

Inhibit  Pulse 

Pulse 

Duration 

Voltage  in 
Normal  State 

1 

+ 

8.k  millisec. 

+0.07 

2 

— 

28 

millisec . 

-0.07 

3 

+ 

1 1 

millisec . 

+0.06 

k 

- 

33 

millisec . 

-0.08 

5 

+ 

9 

millisec . 

+o.ok 

6 

- 

ik 

millisec . 

-0.08 

TABLE  E-1 .  MONOSTABLE  MULTIVIBRATORS  IN  THE 


ITERATIVE  CONTROL  UNIT. 


APPENDIX  F.  SWITCHING  ARRANGEMENT  FOR  INTEGRATORS. 

When  it  is  required  to  switch  the  mode  of  an  integrator, 
one  possible  scheme  is  the  removal  of  the  pot-set  relay  and 
then  placing  a  diode  switch  in  the  summing  junction.  How¬ 
ever,  even  a  small  unbalance  in  the  diode  bridge  after  being 
integrated  will  produce  a  singificant  error.  An  alternate 
scheme  was  devised  to  switch  the  integrator  into  the  initial 
condition,  operate,  and  hold  modes  which  is  not  critically 
dependent  on  diode  bridge  balance.  The  scheme  employed  is 
shown  in  Figure  P-1 . 


I.C 


o 


FIGURE  P-1  .  SWITCHING  ARRANGEMENT  FOR  INTEGRATORS. 

With  the  indicated  applied  voltage,  the  diode  bridge 
will  conduct  heavily  bringing  the  input  voltage  to  the 
integrator  I  very  near  zero  volts.  This  will  correspond  to 
placing  the  integrator  either  in  the  Initial  condition  or  the 
hold  mode.  If  the  diode  switch  has  the  opposite  polarity 
applied,  it  will  be  cut-off  and  I  will  have  the  property 
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of  an  integrator  in  the  operate  mode  with  a  .2 M  input  resistor. 
If  the  applied  polarity  to  the  bridge  is  again  changed, 
causing  it  to  conduct  heavily,  the  integrator  will  now  hold 
its  last  value. 

'The  error  of  I  while  integrating  is  less  than  u.Oli  volts 
in  ten  seconds.  The  drift  while  the  integrator  is  placed 
in  the  hold  mode  is  approximately  .08  volts  in  ten  seconds. 

The  error  however  can  be  minimized  even  further  if  the 
11 gain  of  1U  input  on  integrator  I  is  used,  with  a  gain  of 
1  used  on  integrator  I,  drift  in  the  hold  mode  is  less  than 
0.02  volts  in  1 0  seconds.  It  is  important  to  use  as  large 
a  voltage  as  available  to  put  the  diode  bridge  into  heavy 
conduction.  Tnis  will  help  to  limit  the  drift  when  the 
integrator  is  in  the  hold  mode.  Slightly  better  results  in 
both  the  overate  and  hold  modes  can  be  obtained  if  a 
balanced  diode  bridge  is  used. 


APPENDIX  G.  NEWTON’S  FORWARD  DIFFERENCE  INTERPOLATION 

FORMULA . 

The  Newton' s  Foward  Difference  Interpolation  Formula 
for  ten  points  is: 

fit)  =  a  +  bt  +  ct(t-l)  +  dt ( t-1 ) ( t-2 )  +  et(t-t ) (t-2) (t-3)  + 

2!  3?  IJ1 

gt(t-1)(t-2)(t-3)(t-li)  +  ht  ( t-1  )  (t-2)  ( t  -  3 )  ( t-h ) ( t-9 ) 

5!  Ff 

+  it  (t-1  ) (t-2) (t-3) (t-IQ  (t-S)  ( t-6)  + 

7! 

jt(t-1  )  (t-2)(t-3)(t-)i)(t-5)  (t-6)  (t-7)  + 

kt  ( t-1  )  (t-2)  ( t  —  3  )  (t-lj)  (t-3)  (t-6)  ( t-7  )  ( t— 8 )  ....(G-1) 
9! 

The  values  of  the  stored  function  at  equal  intervals 
of  t  will  be  designated  as  indicated  in  Figure  G-1 . 


FIGURE  G-1 c  NEWTON'S  METHOD  FOR  10  POINTS. 


The  derivatives  of  all  orders  for  expression  (G-1 )  will 


now  be  derived,, 
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f (t)  -  b  +  <2(2t-1)  +  d(3t2-6t+2)  +  e ( 4t3-1 8t^+22t-6 )  + 

2  5  24 

g(5t^-40t3+105t2-10Qt+24)  + 

120 

hd_6t" -75t^+340t3-675t2+54°t-1  20 )  + 

720 

i(7t6-1  26t^+875t^-29l4-0t3+li875t2-3528t+720)  + 

FOI4.O 

ij(8t7-196t6+1Q32t3-9800t’!+27,076t3-39,396t2+26,1  36t-50!^0) 

W7320 

+  k(9tM-288t7+3822t6-27,2l6t^+112,245t1;-269,136t3  +  354,372t^ 
362,880  -21 9, 1 68 t +40, 320) 


f(t)  =  c  +  df  6t-6 )  +  e(12t2-36t+22)  +  g( 20t3-1 20t2+21 0t-1 00 )  + 

F~  2k  120 

h(30t^-300t3+1 020t2-1 350t+548 )  + 

720 

i (42t^«630t^+3500t3-8820t2+9750t-3528 )  + 

^§50 

jJLS^t^-l  1  76t"+9660t^-39  >  200t3+8l  , 228t2-78, 792t+26, 1 36)  + 
40,320 

k(72t7-20l6t6+22,932t-3-1  36,080t^+448,980t  ~807,408t^+ 

362TF8O  708 , ?44t-21 9,168) 


f"(t)  =  d  +  e(2Ut-36)  +  &(60t2-240t+210)  + 

Wi  1 20 

h(120t3-900t2+2,0i4.0t-1350)  + 

720 

i ( 21  0t^-2520t3+1 0,500t2~1 7,640t+9,750)  + 

^0 

.i  (336tr -5880t^+38,640t3-1 1  7, 600t2+1  62,456t-78, 792)  + 

P7320 

k(  50lit^-1  2, 096t^+1 1 4,660 t^-544, 320 t3+346, 940 t2- 

3E2788O  1 ,614,81 6t+709, 744) 
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f(t)  -  e  +  £(_1  20 1- 2)4.0)  +  h(  360t~-1 800t+20k0)  + 

120  720 

i4840t3-7960tP+21  ,000t-1  7.6I4O)  + 

5oKo 

JHM0t''-23,520t3+1 1 9, O20t2-235,200t+1 62,496)  + 
407320 

k(3,02Ut  '-60,480^+498 , 640t^-i  ,632,960t2+ 
362,8730 

2,693 , 880t-1  ,61)4,816) 


=  g  +  h(  720t-1  800 )  +  i(25>20t2-1  9.1  20t+21  . 000  )  + 

720  ^oITo 

l(6720t3-70,660t"+231 ,840t-239,200) 

40,320 

k(l9,120tk-24l ,920t3+1 , 379, 920t2-2, 265, 920t+2, 693,880) 


fVI(t)  =  h  +  i(9040t-1 9,1 20)  4-  .1(20, 160t2-1)i1  ,120t+231  .840)  + 

9040  407320 

k(  60 . 480 1 3 -729. 760 1 ^+2. 781  , 8)40 1- 2, 269, 920 ) 

362,880 


,fv"(t)  =  i  +  j(40,320t-l4l  ,120)  +  k(l8l  ,440t2-1  ,491 ,920t  + 

407320  362,880 

2,791 ,840) 


fvl"  (t)  =  j  +  k( 362, 880t~1  ,491  ,920) 

362,880 


f'*  (t)  =  k 


By  forming  a  table  as  given  in  reference  1 ,  we  may 
find  the  coefficients  a,b,c,...k.  Table  G-1  is  drawn  up 


for  values  a  to  e  in  order  to  demonstrate  the  manner  in 
which  these  coefficients  may  be  found. 
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=a  Ap=b 
Ko 

Kl-K0 

Ki 

K2"K1 

KP 

k3"k2 

k3 

kI|."k3 

k4 


K2“2K1 +ko 
K^-2K2+K>, 

V2K3+K2 


k3-3K2+3Ki-k0 

%-l4.K3+6K2~llK1  +K 

K)r3K3  +  3K2-K1 


TABLE  G-1.  DETERMINATION  OF  COEFFICIENTS  a  to  e, 


Coefficients  a  to  i  are  given  in  Table  G-2  in  terms 
of  the  value  of  the  function  at  the  sampled  points 


a  =  K0 

b  =  *1  -  K0 

c  =  K2  -  210,  +  Kq 

d  =  K3  -  3  Kp  +  3Kn  -  K0 

e  =  Kl  “  I4.K3  +  6K2  -  I4-K1  +  K0 

g  =  K5  -  +  IOK3  "  1uK2  +  ^K1  '  Ko 

h  =  -  6K^  +  15>N[|_  “  2uK3  +  1  5>K2  -  6 10]  + 

i  =  ^  -  7K6  +  21 -  35K^  +  3^K3  -  21 K2  +  7Kn  -  K0 

TABLE  G-2„  COEFFICIENTS  a  TO  i. 

In  order  to  employ  Newton’s  Iterative  Method  the  value 
of  the  derivatives  of  all  orders  of  the  approximation 
function  must  be  known  at  t  =  0.  The  results  for  five  points 
are  given  in  Table  G-3. 
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f  (t) 


a  +  bt  +  _ct  ( t-1  ) 
2! 


O 

+  d(t3-3t~+2t) 
3! 


e  ( t^--6t3+1 1  t2-6t) 

5T 


f(0)  =  a 


+ 


c  (_2 1 - 1  )  +  d[3t2-6t+2)  +  e(l;t3-l8t2+22t-6) 

21  3 !  FT 


f'(0)  =  b  -  .gc  +  ,333d  -  . 2^e 
f(t)  =  c  +  d(6t-6)  +  e(12t2-36t+22) 

it  nr 

f"(0)  =  c  -  d  +  .91 6e 
f(t)  =  d  +  e(2}±t-3& ) 

f (0)  =  d  -  1  ,$e 

i/ 

f(t)  =  e 
f(0)  =  e 


TABLE  G-3.  VALUES  OP  f1(0). 


